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I . INTRODUCTION 


This  document  is  the  Final  Report  provided  by  Systems  Control,  Inc. 

(SCI)  to  the  Office  of  Assistant  Chief  of  Staff  Studies  and  Analysis,  United 
States  Air  Force,  under  contract  F44620-73-C0033,  "Development  of  Methodology 
To  Evaluate  Strategic  Command  and  Control  Systems."  The  project  was  technically 
directed  by  the  Command,  Control  and  Reconnaissance  Division,  Director  of  Strate- 
gic Offensive  and  Defensive  Studies,  Assistant  Chief  of  Staff  Studies  and  Analysis, 
HQ  USAF,  which  provided  valuable  guidance  and  data  throughout  the  conduct 
of  this  research. 

The  principal  objective  of  the  effort  was  to  aid  the  Command,  Control 
and  Reconnaissance  Division  in  the  development  of  methodology  for  evaluating 
the  survivability  and  effectiveness  of  Strategic  Command  and  Control  (C  ) 

Systems.  While  many  issues  must  be  considered  when  analyzing  the  capabilities 
of  Strategic  C2  Systems,  the  principal  one  addressed  in  this  study  involves 
the  communication  systems  which  would  be  employed  in  the  conduct  of  a nuclear 
war.  The  availability  and  time  responsiveness  of  the  communication  systems, 
the  vulnerability  of  the  systems  to  nuclear  and  electronic  warfare  (EW) 
attack,  and  the  relationship  between  communication  system  effectiveness 
and  Strategic  Command  and  Control  effectiveness,  are  examples  of  some  im- 
portant issues  that  the  methodology  addresses. 

The  main  element  of  the  methodology  is  a Strategic  C2  effectiveness 
model  that  evaluates  communication  system  effectiveness  in  terms  of  Strategic 
C2  effectiveness  while  accounting  for  the  salient  factors  which  have  a signifi- 
cant impact  on  system  effectiveness.  The  model  computes  the  statistics  of  message 
arrival  as  a function  of  time  for  a given  network  specification  and  a given 
enemy  nuclear  attack  on  the  network.  The  model  consists  of  two  separate 
but  interrelated  computer  programs  and  it  is  computationally  efficient, 
allowing  Strategic  C2  planners  to  conduct  extensive  parametric  studies 
of  various  offense  tactics  and  defense  communication  systems  and  operating 
procedures.  Factors  that  impact  system  effectiveness  and  that  are  accounted 

for  in  the  model  include : 


Signal  strength  attenuation  resulting  from  atmospheric  nuclear 
detonations 
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Noise  from  enemy  jammers 


• Survivability  (i.e.,  probability  of  being  operable)  of  links  and 
nodes  under  nuclear  attack 


• Time  delays  incurred  by  messages  passing  through  communication 
system  nodes 


Transmission  delays  of  messages  transmitted  on  low  data  rate  links 
and  probability  of  correct  message  receipt  on  these  links 


Directionality  of  communication* links  (i.e.,  one-way  or  two-way 
communication) 


Survivability  of  strategic  forces  at  the  time  of  message  receipt 


An  important  fact  to  note  is  that  communication  system  effectiveness 

11 /i. ’ „ 2 - - . 


(as  well  as  Strategic  C effectiveness)  can  only  be  estimated  in  a statistical 
sense  due  to  the  stochastic  nature  of  the  processes  involved.  The  model 
therefore  mist  be  capable  of  generating  expected  values  and  probability  dis- 
tributions for  all  parameters  of  interest  such  as  the  number  of  nodes  in  a 
specified  class  that  receive  a message  prior  to  a specified  time,  or  the  time 
of  message  reception  at  specified  nodes  within  the  network.  The  probability 
distributions,  which  take  full  account  of  correlated  events  that  inevitably  arise 


in  message  propagation  through  a communication  network,  are  of  particular 
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importance  to  Strategic  C planners . These  distributions  are  necessary  to 

allow  proper  interpretation  of  the  effectiveness  data  generated  by  the  model; 

for  example  consider  the  following  two  probability  distributions:  (1)  all 

20  nodes  of  a class  receive  the  message  by  a specified  time  with  probability 

0.5  or  none  of  the  20  nodes  of  that  class  receive  the  message  by  that  time 

with  probability  0.5;  and  (2)  all  20  nodes  of  the  specified  class  receive  the  message 

by  the  specified  time  independently  with  probability  0.5.  These  two  distribu- 


tions yield  identical  expected  numbers  of  nodes  receiving  the  message  — 10 


yet  would  have  significantly  different  implications  on  Strategic  C2  effectiveness. 
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It  is  also  important  to  note  that  while  a static  connectivity  analysis 

has  historically  been  employed  to  evaluate  communication  system  effectiveness, 

such  an  approach  is  often  misleading  and  inappropriate  for  e raluating  Strategic 

C2  comour ication  system  effectiveness.  The  static  connectivity  approach 

determines,  given  node  and  link  probabilities  of  being  operable,  the  probability 

that  an  operable  path  exists  between  any  two  specified  nodes  at  that  time 

(ignoring  any  link  directionalization  or  delays  along  the  path).  In  a Strategic 

C2  communication  network,  the  important  issue  is  when  a message  arrives  at 

destination  nodes  given  that  it  was  initiated  at  a specified  time  by  a specified 

2 

node.  The  "dynamic  message  propagation"  approach  to  Strategic  C effective- 
ness analysis  described  in  this  report  explicitly  addresses  this  issue  and  can 
be  differentiated  from  the  static  connectivity  approach  as  follows: 

• Link  directionality  is  fully  accounted  for  in  the  dynamic 
approach  whereas  it  is  not  in  the  static  approach. 

• Probability  distributions  of  message  arrival  times  at  specified 
nodes  are  computed  in  the  dynamic  approach  whereas  the  probability 
an  operable  path  exists  between  two  specified  nodes  at  a specified 
time  is  computed  in  the  static  approach. 

• The  relationship  between  the  position  of  a message  in  the  network 
and  the  times  of  link  and  node  destruction  (e.g.,  by  nuclear  attack) 
is  fully  accounted  for  in  the  dynamic  approach  (i.e.,  destruction 

of  a relay  node  after  a message  passes  through  it  does  not  influence 
the  statistics  of  propagation  of  that  message)  whereas  it  is  not  treated 
correctly  in  the  static  approach. 

2 

The  first  point  above  is  important  because  many  links  in  a Strategic  C 

communication  network  are  one-way.  The  second  point  is  important  because, 

2 

for  a variety  of  reasons,  propagation  of  messages  through  a Strategic  C 
communication  network  can  take  a significant  amount  of  time  (i.e.,  minutes 
instead  of  fractions  of  seconds)  and  it  is  end-to-end  message  propagation  time 
that  plays  a critical  role  in  deter&ining  strategic  command  and  control 
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effectiveness.  The  third  point  is  important  because  during  the  time  interval 
a Strategic  C2  communication  system  is  most  needed  and  messages  are  propagat- 
ing through  the  network,  it  is  also  susceptible  to  attack  and  hence  mav  suffer 
outages  at  links  and  nodes;  hence,  when  generating  a candidate  offense  attack 
on  the  network  in  order  to  properly  test  the  network  effectiveness,  one  must 
not  only  determine  which  links  and  nodes  to  attack  but  also  when  the  attacks 
must  of cur  to  be  effective  in  degrading  the  message  arrival  times  at  destina- 
tion nodes. 

Systems  Control,  Inc.  has  contributed  to  the  development  of  the  dynamic 

2 

approach  for  evaluating  Strategic  C effectiveness  by  assisting  the  Command 
Control  and  Reconnaissance  Division  in  formulating  an  appropriate  methodology 
and  by  defining  and  developing  the  two  computer  programs  upon  which  the 
methodology  is  based~the  Network  Status  Model  and  the  Dynamic  Network  Simulator. 
The  Network  Status  Model  computes,  as  a function  of  time,  the  status  of  each 
node  and  link  in  the  network.  Nuclear  effects  calculations,  electronic  warfare 
calculations,  and  hence  all  calculations  relating  to  network  geometry  and 
node /link  vulnerability  are  performed  by  this  model.  The  Dynamic  Network 
Simulator,  which  requires  as  input  the  Network  Status  Model  output,  propagates 
messages  through  the  overall  network  and  computes  the  statistics  relating 
to  Strategic  C2  performance.  Most  of  the  effort  in  this  initial  study  was 
focused  on  the  development  of  the  Dynamic  Network  Simulator;  this  work  will 
be  fully  described  in  succeeding  chapters.  Development  of  the  Network  Status 
Model  is  on-going  and  hence  will  not  be  described  in  great  detail  in  this 
report. 

The  next  chapter  more  fully  discusses  the  problems  of  modelling  and 
analyzing  the  effectiveness  of  Strategic  C communication  networks,  and 
the  methodology  developed  for  addressing  the  problem.  Succeeding  chapters 
describe  the  methodology  in  more  detail. 
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II . THE  STRATEGIC  C2  PROBLEM  AND  SOLUTION  METHODOLOGY 

The  Strategic  C2  networks  of  concern  here  link  the  National  Command 
Authority  (NCA)  to  the  strategic  forces — at  present  consisting  of  the  Strategic 
TRIAD:  Minuteman  Missiles,  SAC  Bombers,  and  Polaris  Submarines.  The  Strategic 

C2  communication  problem  is  concerned  with  the  development  and  operation  of 
communication  systems  to  be  employed  in  the  conduct  of  a general  nuclear  war. 

The  problem  may  be  broken  into  three  subproblems:  the  attack  warning  problem-- 

determination  that  an  enemy  attack  has  been  initiated;  the  delivery  of  attack 
orders  problem — transmission  of  the  Emergency  Action  Message  (EAM)  from  the  NCA 
to  the  strategic  forces;  and  the  battle  management  problem — the  exchange  of 
messages  between  the  NCA  and  the  Strategic  forces  during  the  course  of  the  war. 
This  report  is  primarily  concerned  with  modelling  and  evaluation  of  communica- 
tion systems  employed  in  delivering  attack  orders. 

2 

Description  of  the  Strategic  C Communication  Network 

A Strategic  C network  can  consist  of  a variety  of  land-based,  airborne, 
seaborne  and  satellite  nodes  electromagnetically  linked.  Potentially,  the  net- 
work may  nave  land  lines,  or  electromagnetic  links  in  any  of  the  bands  ELF-LF, 

MF-HF,  and  VHF-EHF.  However,  the  low  frequency  links  (ELF-LF)  and  the  high 
frequency  links  (VHF-EHF)  are  most  likely  to  survive  a nuclear  attack  and 
hence  only  these  links  will  be  explicitly  considered  here.  However,  the  method- 
ology developed  herein  is  also  applicable  to  the  other  forms  of  links.  Because 
the  nodes  and  links  of  a Strategic  C network  are  subject  to  possible 
destruction  and  disruption  by  nuclear  bursts  and  enemy  jamming,  a multiplicity 
of  paths  are  generally  provided  between  the  NCA  and  the  strategic  forces  in 
particular,  two  nodes  are  often  connected  by  two  or  more  links  of  differing 
types. 

Properties  of  Nodes  and  Links.  The  propagation  of  the  EAM  through  a Strategic 
C2  Network  depends  upon:  (1)  the  survivability  of  nodes  and  their  times  of 
destruction,  (2)  the  time  delay  at  each  node  between  receipt  and  retransmission 
of  the  message  due  to  nodal  processing  procedures  such  as  decoding  and  coding. 
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and  (3)  the  delay  for  each  link,  that  is,  the  time  from  initial  transmission 
of  the  message  at  the  sending  node  until  it  is  completely  and  correctly 
received  at  the  receiving  node.  Link  delays  are  a function  of  the  length  of 
the  EAM,  the  form  of  link,  and  the  probability  of  incorrect  reception  of  the 
message.  For  example,  the  transmission  time  at  ELF- LF  may  be  minutes  and 
a message  may  be  transmitted  several  times  before  it  is  correctly  received; 
hence  the  link  delay  may  be  a multiple  of  the  message  transmission  time  or 
even  infinity  if  the  link  is  broken.  The  message  transmission  times  at  VHF-EHF 
are  much  smaller  but  still  not  negligible. 

Types  of  Nodes.  A node  may  be  moving  or  stationary  and  may  be  classified  as 

being  in  one  of  three  classes:  source  nodes,  ordinary  nodes,  and  special 

nodes.  A source  node  is  a node  at  which  the  EAM  is  initiated;  typically 

2 

there  will  be  only  one  source  node  in  a Strategic  C network.  An  ordinary 
node  is  neither  the  source,  nor  a destination,  of  the  EAM;  it  serves  to  relay 
the  message  and  it  is  always  linked  to  other  nodes  from  which  it  can  receive 
the  message  and  to  which  it  can  send  the  message.  A special  node  is  any  node 
for  which  the  user  wants  statistics  on  arrival  of  the  EAM. 

Types  of  Links.  A link  between  two  nodes  may  be  one-way  or  two-way,  and  moving 
or  stationary  depending  on  whether  or  not  one  of  the  two  nodes  it  connects 
is  moving.  For  present  purposes  it  is  convenient  to  separate  all  links  into 
four  classes:  ELF-LF  links,  VHF-EHF  links,  satellite  links,  and  special  links. 

ELF-LF  and  VHF-EHF  links  are  stationary  (aircraft  and  seaborne  nodes  are 
assumed  stationary  for  time  and  distance  scales  of  interest  here),  while 
satellite  and  special  links  are  moving  VHF-EHF  links.  The  first  three  types 
of  links  are  assumed  to  be  operating  prior  to  injection  of  the  EAM  into  the 
network.  Special  links  correspond  to  nodes  which  move  a great  distance  during 
the  simulation  and  do  not  become  operable  until  a fixed  time  after  the  sending 
node  of  the  link  receives  the  message.  Special  links  may  be  treated  in  a 
manner  similar  to  fixed  VHF-EHF  links  by  creating  fictitious  links  and  nodes: 
the  message  travels  instantaneously  over  a fictitious  perfect  link  to  a 
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fictitious  node  which  than  becomes  a sending  node  for  the  special  link  after 
a fixed  delay  corresponding  to  the  time  required  for  the  link  to  become 

operable. 

Decomposition  of  the  Problem. 

The  status  of  the  nodes  (i.e.,  dead  or  alive)  and  the  status  of  the  links 
(i.e.,  their  ability  to  transmit  the  EAM)  at  any  given  time  is  a function  of 
certain  external  events  (such  as  nuclear  bursts  ana  jamming)  caused  by  complex 
physical  phenomena,  while  the  progress  through  the  network  of  the  EAM  given  node 
and  link  status  is  statistical.  Therefore  it  is  convenient  and  desirable  to  de- 
compose the  problem  of  evaluating  Strategic  C2  effectiveness  into  two  parts: 
determination  of  network  status  as  a function  of  time  followed  by  determination 
of  network  effectiveness  in  propagating  the  EAM  to  the  strategic  forces.  This 
decomposition  is  possible  because  while  network  effectiveness  depends  upon  the 
time  history  of  network  status,  the  opposite  is  not  true.  The  separation  between 
network  status  and  network  effectiveness  calculations  is  somewhat  arbitrary  and 
should  be  based  on  convenience  and — since  the  calculations  will  be  made  on  a 
digital  computer — on  computer  efficiency  considerations. 

Determination  of  Network  Status.  To  determine  the  network  status  one  must 
compute  the  status  of  each  node  and  each  link  at  each  time  increment  of  interest. 
(The  assumption  is  made  here  that  time  may  be  quantized  into  fixed  increments.) 

An  important  question  is  whether  one  should  compute  the  status  of  all  nodes 
and  links  at  a given  time  and  cycle  through  the  times  of  interest  or  whether 
one  should  compute  the  status  of  a given  node  or  link  for  all  times  and  cycle 
through  the  nodes  and  linkB.  Since  the  computations  made  to  determine  node 
or  link  status  at  a given  time  may  provide  information  that  is  useful  to  determine 
the  status  of  that  node  or  link  at  a later  time,  the  latter  organization  is 
preferable.  Figure  2.1  is  a high-level  flow  chart  of  a computer  program  for 
determining  network  status  (the  Network  Status  Model)  utilizing  this  organization. 
Table  2.1  is  a list  of  inputs  required  by  this  model,  and  Table  2.2  is  a list  of 
outputs  from  this  model.  Link  and  node  status  definitions  and  calculations  will 
be  more  fully  discussed  in  later  sections. 
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TABLE  2.1  INPUTS  TO  THE  NETWORK  STATUS  MODEL 


I 


I 


i 


i 


I 


\ 


1.  Node  Information 

a.  Location:  latitude,  longitude,  altitude  of  each  node. 

b.  Hardness:  number  for  fixed  nodes,  (time,  cumulative  probability 

of  destruction)  pairs  for  other  nodes. 

c.  Node  Delay:  Minimum  node  delay  and  Interval  of  node  delays  for  all 

nodes.* 

d.  Node  Class  for  each  Special  Node.* 

e.  Receipt  Value  vs.  Time:  (time,  value  of  message)  pairs  for  each 
special  node.* 

2.  Link  Information 

a.  Network  Geometry:  sending  node  and  receiving  node  for  all  links. 

b.  Link  Class:  either  ELF-LF,  VHF-EHF,  satellite  or  special  for  each 

link. 

c.  Message  Transmission  time  for  all  links.* 

d.  Time  to  Start  of  Message  Transmission  for  all  special  links.* 

e.  Parameters:  transmitter  effective  radiating  power,  midband  frequency, 

transmitting  and  spread  bandwldths,  modem  type,  data  rate,  receiving 
antenna  gain  and  temperature,  receiver  line  loss,  side  lobe  gain, 
main  lobe  beamwldth,  number  of  characters  In  digital  alphabet,  and 
number  of  phases  for  each  llnk.^ 

3.  Burst  Information 

a.  Location:  latitude,  longitude,  altitude  of  each  burst. 

b.  Time  of  each  burst. 

c.  Parameters:  yield,  fission  fraction,  and  black  body  radiating 

temperature  for  each  burst,  CEP  of  the  warhead. 


* These  inputs  are  not  used  by  the  Network  Status  Model  but  merely  transferred 
to  the  outputs  for  use  by  the  Dynamic  Network  Simulator. 

t This  is  a preliminary  list  of  link  parameters.  Because  of  the  large  number 
of  parameters  and  links  a small  number  of  link  types  will  be  defined 
and  parameters  stored  in  the  model  for  each  type;  input  will  merely  specify 
the  link  type. 
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a.  Location:  latitude,  longitude,  altitude  of  each  jammer. 

b.  Parameters:  effective  radiated  power,  midband  frequency,  and 

bandwidth  for  ea^n  jammer. 


TABLE  2.2  OUTPUTS  OF  THE  NETWORK  STATUS  MODEL 
(Inputs  Required  by  the  Dynamic  Network  Simulator) 

1.  Node  Information 

a.  Node  Status  vs.  Time:  (time,  cumulative  probability  of  destruction) 

pairs  for  each  node. 

b.  Node  Delay:  Minimum  node  delay  and  interval  of  node  delay  for  all 

nodes.* 

c.  Node  Class  for  each  special  node.* 

d.  Receipt  Value  vs.  Time:  (time,  value  of  message)  pairs  for  each 

special  node.* 

2.  Link  Information 

a.  Network  Geometry:  sending  node  and  receiving  node  for  each  link.* 

b.  Link  Class:  either  ELF-LF,  VHF-EHF  satellite,  or  special  for  each 

link.* 

c.  Message  Length  for  each  ELF-LF  link.* 

d.  Time  to  start  of  Message  Transmission  for  each  special  link.* 

e.  Link  Status  vs.  Time:  (time,  probability  of  acceptable  message)  pairs 
for  each  ELF-LF  link,  (time,  probability  of  availability)  pairs 

for  each  other  link. 


* These  outputs  of  the  Network  Status  Model  are  inputs  to  it  and  are  unchanged 
by  it. 
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Determination  of  Network  Effectiveness.  The  propagation  of  an  EAM  through  a 
Strategic  C2  network  is  stochastic  because  node  destruction,  node  delays, 
and  transmission  delays  are  random  processes.  (A  description  of  the  models 
developed  for  these  random  processes  is  given  in  Section  III.)  Network  effectiveness 
is  determined  from  the  statistics  of  message  arrival  at  the  special  nodes. 

(A  discussion  of  appropriate  statistics  is  given  in  Section  V.)  Two  general 

methods  exist  for  determining  the  statistics  of  a stochastic  process:  statistical 

simulation  and  Monte  Carlo  simulation.  In  a statistical  simulation  the 

probability  distribution  of  the  state  of  the  process  is  propagated  forward 

in  time;  whereas  in  a Monte  Carlo  simulation  a number  of  runs  are  made  with 

the  random  processes  simulated  by  appropriate  use  of  random  number  generators. 

Unfortunately,  for  all  but  the  simplest  processes  the  number  of  states  of  the 

system— hence  the  dimension  of  the  probability  distribution— can  be  very 

large  rendering  stochastic  simulation  impractical.*  Fortunately,  use  of 

the  shortest  path  algorithm  (as  detailed  in  Section  IV)  permits  each  run 

of  a Monte  Carlo  simulation  to  be  made  efficiently  so  that  sufficient  runs 

may  be  made  to  obtain  reliable  statistics  with  a reasonable  amount  of  computer 

time  (the  number  of  runs  required  is  discussed  in  Section  V).  Figure  2.2 

is  a high  level  flow  chart  of  a computer  program  (the  Dynamic  Network  Simulator) 

for  determining  network  effectiveness  utilizing  Monte  Carlo  simulation. 

Table  2.2  is  a list  of  inputs  required  by  this  model,  and  Table  2.3  is  a list 
of  outputs  generated  by  this  model.  For  each  node  a determination  of  node 
destruction  time  (possibly  infinity)  is  made  prior  to  message  propagation. 

On  the  other  hand,  to  avoid  computing  node  and  link  delays  at  every  time 
increment,  when  in  fact  they  are  needed  only  at  the  time  the  message  arrives 
at  the  node  or  link,  these  delays  are  computed  within  the  message  propagation 
algorithm. 

Efficient  Use  of  Network  Status  Model  and  Dynamic  Network  Simulator.  Because 
the  Network  Status  Model  contains  nuclear  effects  codes  and  generally  requires 


A 

* For  each  node  of  a strategic  C network,  at  a given  time  it  must  be  known  if 
and  when  a message  has  arrived.  Thus,  at  time  k each  node  can  be  in^any  one 
of  k + 1 states  and  a system  of  m nodes  can  be  in  any  one  of  (k  + 1)  states— 
an  astronomical  number  for  typical  values  of  k and  m. 
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Propagate  message  through  network 
determining  node  delay  (from  node  clase 
and  random  number)  and  link  delay  for 
each  link  leaving  the  node  (from  link 
atatua  va  time,  link  claaa,  and  random 
number)  at  time  message  arrlvea  at  node 


Update  desired 
atatlatlcs 


FIGURE  2.2  DYNAMIC  NETWORK  SIMULATOR 
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TABLE  2.3  OUTPUTS  OF  THE  DYNAMIC  NETWORK  SIMULATOR* 

1.  Time  Period  Histograms  (See  Section  V)  . 

2.  Arrival  Tim*  Distributions  (See  Section  V). 

3.  Probability  of  any  Event  that  can  be  defined  as  occurring  or  not 
occurring  on  a Monte  Carlo  Run. 

4.  Probability  Distribution  Over  a Set  of  Mutually  Exclusive  Events 
(Outputs  1,  2,  and  3 are  special  cases  of  output  4). 

3.  Expected  Value  of  Performance  Index  (See  Section  V). 

6.  Expected  Value  of  any  Random  Variable  whose  value  may  be  defined  for 
each  Monte  Carlo  Run  (Outputs  3 and  5 are  special  cases  of  output  6). 

7.  For  any  Monte  Carlo  run  a map  of  the  dominant  Nodes  and  Links  that 
actually  carry  the  message. 

8.  For  any  link,  the  probability  it  will  carry  the  message  (a  special  case 
of  output  3) . 

9.  For  any  node,  the  probability  it  will  relay  the  message  (a  special  case 
of  output  3). 


* This  list  is  not  meant  to  exhaust  the  possibilities  of  outputs,  but  rather 
to  suggest  some  of  them.  Any  output  summarizing  the  results  of  all  Monte 
Carlo  runs  may  be  viewed  as  a set  of  outputs  of  Type  6. 
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a significant  amount  of  computer  time,  while  the  Dynamic  Network  Simulator 
is  computationally  much  faster,  it  is  important  to  utilize  the  two  models 
in  an  efficient  manner.  Typically,  one  will  be  interested  in  a baseline 
Strategic  network  and  several  proposed  modifications  of  the  baseline.  To 
perform  trade-off  studies  efficiently  it  is  desirable  to  develop  an  augmented 
baseline  network  such  that  the  baseline  and  each  of  the  modified  baselines  is 
a subnetwork  of  the  augmented  baseline  network.  The  Network  Status  Model  is 
then  used  to  obtain  the  status  of  the  augmented  baseline  network  and  a series 
of  simulations  are  made  with  the  Dynamic  Network  Simulator  in  which,  for  each 
simulation,  only  those  nodes  and  links  comprising  the  system  of  current  interest 
are  considered.  Furthermore,  it  may  be  desirable  to  perform  a large  number  of 
trade-off  studies  with  the  Network  Status  Model  embodying  greatly  simplified 
nuclear  codes  followed  by  selective  use  of  detailed  nuclear  codes  to  verify 
key  results  obtained  with  the  simplified  codes. 


III.  DETERMINATION  OF  NODE  AND  LINK  PROPERTIES 


In  the  previous  section  it  was  pointed  out  that  propagation  of  the  EAM 
2 

through  the  Strategic  C network  is  dependent  upon  how  long  each  node  sur- 
vives, the  message  delay  at  each  node,  and  the  transmission  delays  along  each 
link.  This  section  describes  how  these  quantities  are  determined. 

Node  Properties 

Two  node  properties  are  of  interest:  node  survival  and  node  delay. 

These  properties  are  individually  discussed  below. 

Node  Survival.  Consider  a given  node.  For  each  nuclear  burst  in  the  vicinity 
of  the  node  there  is  a non-zero  probability  the  node  will  be  made  inoperable. 
This  probability  depends  upon  the  size  and  location  of  the  burst  and  the 
hardness  and  location  of  the  node.  A model  for  determining  this  probability 
(for  over-pressure  kill  only)  is  given  in  Appendix  A.  By  use  of  this  model 
a number  of  (time,  probability  of  destruction)  pairs  (t±,  q±)  may  be  generated 
with  each  time  t±  corresponding  to  each  burst  that  has  a non-negligible 
probability  qi  of  destroying  the  node  at  time  t^  ^ < t2  < From 

the  pairs  (t^,  q^)  a set  of  time-cumulative  probability  of  destruction  pairs 
(tjt  q±)  can  be  generated—where  q*  is  the  probability  that  the  node  is 
destroyed  by  the  burst  at  time  ^ or  by  any  other  burst  prior  to  time  t . 
The  q^  and  q*  obey  the  following  relationships: 

qi*  - V 

= + 0-qt*)  qi+1  for  i >_  1 (3-1) 

To  determine  the  time  of  destruction  of  a node  a random  unit  (i.e.,  random 
number  uniform  on  [0,1])  is  generated  and  compared  with  the  q*  s.  The  time  of 
destruction  is  the  earliest  t^  for  which  q*  is  greater  than  the  random  number. 
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Computation  of  q^  the  probability  of  destruction  by  a single  burst, 

is  performed  in  the  Network  Status  Model.  The  calculations  represented  by 

(3-1)  can  be  made  logically  in  either  the  Network  Status  Model  or  in  the 

Dynamic  Network  Simulator — a somewhat  arbitrary  decision  was  made  to  perform 

these  calculations  in  the  Network  Status  Model.  In  the  present  version  of 

the  Network  Status  Model  (t,,  q*)  pairs  are  computed  only  for  ground-based 

11  * 

nodes — the  user  must  specify  (t  , q*)  pairs  for  other  node  types  (e.g.,  satellite 

1 1 2 

and  aircraft  nodes) . During  each  Monte  Carlo  simulation  of  the  C system 
the  time  of  node  destruction  is  determined  by  the  Dynamic  Network  Simulator 
prior  to  use  of  the  shortest  path  algorithm. 


Node  Delay 


mum 


Node  delays  are  currently  input  for  each  node.  They  are  assumed  to  be 
uniformly  random  between  minimum  and  maximum  values  ^t^  and 

(*  ^^inimum  + ^Interval)*  actua*  noc^e  delay  At  is  determined  for 

each  Monte  Carlo  run  in  the  Dynamic  Network  Simulator  by  selecting  a random 
unit  r for  each  node  prior  to  the  time  it  is  to  relay  the  EAM.  At  the 
time  it  is  to  relay  the  message  the  delay  is  determined  from 


At 


^Minimum  + r ^Interval  ^ 


(3-2) 


For  each  fictitious  node  generated  by  a special  link,  the  node  delay  is  equal 
to  the  time  after  message  receipt  that  the  link  becomes  active  (see  Section  II)  . 


Link  Properties 

It  is  assumed  that  the  received  signal  strength  and  noise  can  be  modelled 
by  a Gaussian  distribution.  Accordingly,  the  moan  and  standard  deviation  of 
signal  and  noise  at  the  receiving  nodes  must  be  computed.  For  ELF-LF  and 
VHF-EHF  links  these  values  are  used  in  different  manners  to  compute  link 
delays . 

Signal  and  Noise.  The  signal  strength  at  the  receiving  node  of  a link  in 
the  absence  of  nuclear  effects  is  a function  of  sending  and  receiving  node 


* There  are  standard  procedures  cor  doing  this. 
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locations  and  the  link  class  (which  determines  transmitter  and  receiver 
characteristics,  frequency,  etc.).  Nuclear  bursts  create  absorbing  regions 
(the  only  effect  thus  far  considered)  that  lower  the  signal  strength  at  the 
receiver.  The  amount  of  signal  attenuation  is  a function  of  burst  size  and 
location  m addition  to  the  link  iocation  and  class.  Noise  at  the  receiver 
is  the  sum  of  ambient  (atmospheric)  noise  (a  function  of  link  class)  and 
noise  from  jammers  (a  function  of  jammer  location  and  type,  and  link  location 
and  class).  Description  of  the  models  required  to  determine  nuclear-free 
signal  strength,  nuclear  attenuation  levels,  and  ambient  and  janming  noise 
levels  is  beyond  the  scope  of  this  report.  Suffice  it  to  say  that  the  requisite 
computations  are  to  be  made  in  the  Network  Status  Model. 

ELF-LF  Links.  Given  a particular  value  of  the  Signal  to  Noise  ratio  (S/N) 
there  is  a resulting  Probability  of  Character  Error  (PCE)  which  is  a charac- 
teristic of  the  type  of  communication  link.  Given  a particular  value  of  PCE 
there  is  a resulting  probability  that  the  message  will  be  accepted  (PAM) 
which  is  a function  of  how  the  message  is  structured  into  characters  and  the 
error  correction  procedures  that  are  employed.  Thus  there  is  a one  to  one 
relation  between  S/N  and  PAM  which  is  defined  by  series  application  of  the 
two  curves  relating  PCE  to  S/N  and  PAM  to  PCE.  If  the  S/N  varies  over  the 
message  transmission  interval,  it  is  appropriate  to  use  the  smallest  S/N  value 
over  the  interval  since  even  a short  drop  in  S/N  may  result  in  many  character 
errors  which  render  the  overall  message  unacceptable.  Thus  starting  from  the 
Gaussian  distribution  for  S/N  it  is  straightforward  to  obtain  a curve  nf 
PAM  versus  time,  which  yields  the  probability  that  a single  transmission  of 

the  message  starting  at  a given  time  will  be  acceptable  when  it  is  completely 
received. 

Consider  a given  ELF-LF  link  and  assume  that  a single  transmission  of  the 
EAM  takes  At  time  units,  that  the  EAM  is  first  transmitted  at  time  t,  and 
that  it  is  sent  repeatedly  until  it  is  correctly  received  at  the  receiving 
node  or  the  final  time  of  interest,  tf,  has  passed.  Clearly,  the  message 

may  be  correctly  received  at  t + At,  t + 2At t + (n-1)  At,  t + nAt  - t : 

therefore  the  possible  link  delays  are  At,  2At nAt,  where  a delay  of  f 

(n+1)  At  implies  that  the  message  never  reaches  the  receiving  node  (during 
the  time  interval  of  interest).  To  determine  the  link  delay  which  actually 
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applies  on  a given  Monte  Carlo  run,  the  cumulative  probability  distribution 
p*,...,  p*  over  the  delay  values  can  be  determined  from  the  curve  of  PAM 
versus  tiHe  (using  an  equation  similar  to  3-1)  and  compared  with  a random 

unit. 


Determination  of  the  cumulative  probabilities,  p*,  may  be  made  more 
complex  if  certain  error  correcting  procedures  are  employed.  For  example, 
in  one  procedure  the  message  could  be  correctly  received  at  a given  time  on 
the  basis  of  the  last  transmission,  the  last  three  transmissions,  or  the  last 
five  transmissions.  The  ability  to  correctly  receive  a message  is  then  defined 
for  each  link  class  to  be  a set  of  three  curves,  one  for  each  of  1,  3,  and  5 
transmissions,  relating  the  probability  of  correctly  receiving  the  message 
to  PCE.  Since  PCE  varies  with  S/N  which  in  turn  varies  with  time,  the  minimum 
value  of  S/N  over  the  interval  of  interest  should  be  used.  For  example:  the 


probability  pk(l)  that  the  message  was  correctly  received  at  time  t + kAt 
on  the  basis  of  the  kth  transmission  alone  is  found  by  determining  PCE  using 


S/N 


minimum 


S/N 


(3-3) 


i such  that  t + (k-DAtft^t  + kAt 


and  setting  pk(l)  equal  to  the  probability  of  correct  message  receipt  cor- 
responding to  this  PCE  using  the  single  transmission  curve.  The  probabilities 
Pk(3)  and  pfc(5)  that  the  EAM  was  correctly  received  at  time  t + kAt  on  the 
basis  of  transmissions  k-2,  k-1,  and  k or  k-4,  k-3,  k-2,  k-1,  and  k, 
respectively,  are  determined  similarly  except  that  the  minimizations  of  (3-3) 
should  be  carried  out  over  the  intervals  t + (k-3)At  or  t + (k-5)At  to  t + kAt 
and  the  appropriate  curves  used.  An  equation  similar  to  (3-1)  can  then  be 
used  to  determine  the  cumulative  probabilities  p*  from  pk(l),  Pk(3) , and  P^(5)‘ 
Figure  3-1  is  a flow  chart  of  the  complete  procedure  used  to  determine  the 

p*'s  from  the  (time,  S/N  pairs, 
k i 


If  the  data  transferred  between  the  Network  Status  Model  and  the  Dynamic 
Network  Simulator  were  the  time-varying  S/N,  the  logic  of  Figure  3-1  could  be 
carried  out  as  part  of  the  shortest  path  algorithm  as  follows.  A random  unit 


could  be  selected  for  each  link  prior  to  determination  of  the  p*'s  at  the  time 


Determine  y, (J)  from 
(t.,  S/N^)  pairs  using 
appropriate  curves 


at  which  the  message  is  to  be  sent  over  the  link  • Each  time  a new  p*  is 
computed  it  would  be  compared  with  the  random  unit  and  when  p*  exceeded  the  ran- 
dom unit,  the  message  would  be  declared  to  be  correctly  transmitted  with  a link 
delay  kAt.  A modified  procedure  is  needed  if  the  data  transfer  between  the  NSM 
and  the  DNS  is  the  time-varying  values  of  probability  of  acceptable  message 
(PAM), 


VHF-EHF  Links . Because  the  EAM  requires  only  a short  transmission  time  at 
VHF-EHF,  the  links  are  modelled  in  a simpler  fashion  than  the  ELF-LF  links 
From  the  signal  and  noise  strength  time  histories,  a time  history  of  the  pro- 
bability P of  successful  transmission  is  determined  in  the  Network  Status 

Model  in  the  form  of  (time,  P ) pairs  (t . , P ).  (Again,  the  means  of  per- 

t i t^ 

forming  this  determination  is  beyond  the  scope  of  this  report.)  These  pairs 
are  utilized  in  the  Dynamic  Network  Simulator  as  follows : For  each  VHF-EHF 

link  a random  unit  is  generated  prior  to  the  time  the  EAM  is  to  be  sent  over 
the  link.  At  the  time  ti  that  the  message  is  sent  over  the  link,  the  random 
unit  is  compared  with  Pfc  corresponding  to  the  largest  ^ less  than  t.  The 
message  is  successfully  deceived  depending  on  whether  or  not  the  random  unit 
exceeds  P.  . 


IV.  SIMULATION  OF  MESSAGE  PROPAGATION 


Given  the  statistics  describing  the  node  destruction  times,  the  node 

delays,  and  the  link  delays,  the  propagation  of  the  EAM  throughout  the 
2 

Strategic  C Network  must  next  be  modelled.  Section  II  Indicated  that  this 
is  accomplished  through  use  of  the  shortest  path  algorithm  (SPA).  A descrip- 
tion of  this  procedure  is  the  major  subject  of  this  section,  but  first  the 

I 

motivation  for  using  the  SPA  is  discussed  briefly. 


Possible  Methods  of  Analysis 

Connectivity  Algorithms.  The  classical  tools  for  analyzing  networks  subject 

to  node  and  link  failures  are  connectivity  algorithms  which  proceed  as 

follows:  based  on  probabilities  of  failure  and  random  units,  each  node  and 

link  is  determined  to  be  either  operable  or  Inoperable.  A connectivity  algorithm 

then  efficiently  partitions  the  network  into  subsets  of  connected  nodes  such 

that  communication  between  node  A and  node  B is  possible  if  A and  B are  in 

the  same  subset.  By  performing  a Monte  Carlo  sequence  of  runs,  statistics 

such  as  the  probability  that  A can  communicate  with  B,  or  other  statistics 

relating  to  the  connectivity  of  nodes  in  the  network,  can  be  determined. 

Reference  1 gives  a detailed  discussion  of  this  approach  as  it  applies  to 

2 

analysis  of  Strategic  C Networks.  Unfortunately,  connectivity  algorithms 
suffer  two  shortcomings  when  applied  to  the  problem  of  delivering  attack 
orders.  First,  they  only  apply  to  two-way  links — if  A is  connected  to  B then 
B is  connected  to  A;  secondly,  they  Ignore  dynamics — the  times  at  which  links 
or  nodes  become  inoperable  are  not  taken  into  account  and  the  node  and  link 
delays  are  limited  to  the  values  of  zero  and  infinity.  Connectivity  algorithms 
are  more  appropriate  for  steady  state  analyses  of  a network  in  which  many 
messages  are  being  sent  and  in  which  elements  become  operable  and  inoperable 
at  random  times.  Attempts  to  extend  the  connectivity  ideas  to  treat  the 
delivery  of  attack  orders  problem  proved  fruitless  because  of  these  difficulties. 

2 

Shortest  Path  Algorithm.  In  studying  the  Strategic  C mission  of  delivering 
attack  orders  the  basic  assumption  is  made  that  each  node,  after  an  appropriate 
delay  following  message  receipt,  transmits  the  EAM  along  all  links  leaving  it. 


22 


With  this  procedure  it  is  clear  that  the  message  will  arrive  at  every  node 

in  the  network  at  the  earliest  possible  time.  The  shortest  path  algorithm 
(SPA) , known  to  network  theorists  (see  Ref enence  2) , can  be  applied  to  solve 
for  the  resulting  message  arrival  times.  The  path  length  of  each  link  is 
defined  to  be  the  sum  of  the  link  transmission  delay  and  the  delay  at  the 
sending  node.  Since  the  SPA  is  a very  efficient  algorithm,  it  is  quite 
suitable  for  ncdelling  GAM  propagation  in  a Strategic  C network.  It  should 
also  be  noted  that  while  the  SPA  can  be  used  to  model  dynamic  message  propaga- 
tion it  can  also  be  used  to  determine  static  connectivity.  This  is  done  simply 
by  setting  delays  either  to  zero  or  to  infinity;  see  Appendix  B for  details. 

Use  of  the  Shortest  Path  Algorithm 

A detailed  description  of  the  SPA  as  applied  to  the  delivery  of  attack 
orders  problem  will  be  presented,  followed  by  an  example  of  its  use.  For 
the  original  derivation  of  the  SPA  the  reader  is  referred  to  Reference  2. 

Description  of  the  SPA.  The  SPA  is  a procedure  for  labeling  each  node  with 
the  length  of  the  shortest  path  between  it  and  the  source  node.  (If  path 
length  is  time,  the  label  for  each  node  is  the  time  it  receives  the  EAM.) 
Suppose  some  of  the  nodes  have  been  so  labeled  and  refer  to  these  nodes  as 
"considered".  Further  suppose  that  each  of  the  "unconsidered"  nodes  is 
labeled  with  the  earliest  time  that  it  can  receive  the  message  over  a single 
link  from  a "considered"  node.  Now  take  the  "unconsidered"  node  whose 
label  is  smallest.  Clearly  it  cannot  receive  the  message  any  sooner  than  its 
label  from  any  of  the  other  "unconsidered"  nodes;  therefore  it  is  appropriately 
labeled  and  may  be  removed  from  further  consideration  except  that  the  label 
of  each  remaining  "unconsidered"  node  must  be  updated  by  examining  each  link 
leaving  it.  This  procedure  of  picking  the  "unconsidered"  node  whose  label  is 
smallest  is  initiated  with  all  nodes  "unconsidered",  the  source  node  labeled 
with  the  EAM  injection  time,  and  the  remaining  nodes  labeled  infinity.  The 
procedure  is  repeated  until  all  nodes  have  been  "considered". 
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Figure  4 .1  is  a flow  chart  of  the  SPA.  The  source  node  is  N^  ^ is 
the  ith  node,  tj  is  the  time  of  destruction  for  N^  Lj  is  the  j link,  d^  is 
the  node  delay  for  plus  the  link  delay  for  (computed  in  the  SPA  as 
described  in  Section  III),  and  Nll(j)  is  the  node  receiving  the  EAM  from  ^ 
over  L . For  convenience  two  lists  are  maintained  by  the  program:  a consid- 

eration list,  in  which  all  "unconsidered"  nodes  with  finite  label  are  listed 
in  order  of  increasing  label,  and  a node  list  in  which  the  label  of  each  node 
is  listed.  Note  that  in  Box  10  a check  is  made  to  insure  that  both  sending 
and  receiving  nodes  will  be  operable  until  message  transmission  is  complete. 

(Note  further  that,  if  k*(i')  is  finite,  only  the  sending  node  need  be  examined 
since  tj,  > k*(i')  and  only  k < k*(i)  are  of  interest.) 

A Simple  Example.  To  illustrate  the  use  of  the  SPA  an  example  with  two  Monte 
Carlo  rims  will  be  presented.  The  network  under  consideration  in  this  example 
is  given  in  Fig.  4.2.  The  goal  is  to  transmit  a message  from  Nj  to  N3» 

Nominally  and  are  fast  links  with  a delay  of  1 and  a slow  link  with 
a delay  of  3.  The  message  may  not  be  received  correctly  at  the  first 
transmission  from  a node;  therefore  on  each  link  from  the  node  the  message 
is  sent  over  and  over  until  it  is  correctly  received  by  the  receiving  node.  Given 
the  probability  of  correct  reception  for  eacn  transmission  and  the  nominal  delay 
for  each  link,  the  distribution  of  the  delays  along  the  link  is  readily 
determined.  For  this  example  these  delays  are  taken  to  be  independent.  At 
time  3 a burst  occurs  that  kills  node  2 with  probability  y. 


FIG.  4.2  EXAMPLE  NETWORK 
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For  Monte  Carlo  Run  1,  a coin  was  flipped  to  see  if  N2  was  destroyed  at 
time  3 and  it  was  found  that  the  node  was  not  destroyed.  The  iterations 
taken  by  the  shortest  path  algorithm  are  given  in  Table  4.1.  During  iteration 
0,  the  algorithm  was  initialized  by  setting  k*(l)-l,  k*(3)  ■ 00  and  placing 
on  the  consideration  list  at  level  1.  At  the  start  of  iteration  1,  the  top 
node  on  the  consideration  list  was  N^.  This  node  was  deleted  and  the  delays 
along  and  were  computed  using  the  appropriate  distributions  as  given 
in  Table  4.2.  Two  random  numbers  were  selected;  both  were  less  than  0.667 — 
therefore,  d^  was  set  equal  to  1,  d^  was  set  equal  to  3,  N2  was  placed  on 
the  consideration  list  at  level  2,  and  at  level  4. 

At  the  beginning  of  iteration  2,  was  at  the  top  of  the  consideration 

list;  it  was  deleted  and  d23  selected  via  Table  4.3.  A random  number  was 
selected  that  lay  between  0.889  and  0.907;  therefore  d23  was  set  equal  to  3. 
Since 


k*(2)  +d23«2+3>4«  k*<3), 

the  position  of  on  the  consideration  list  and  the  value  of  k*(3)  were 
left  unchanged.  At  the  start  of  iteration  3,  only  N3  is  on  the  consideration 
list.  Because  it  has  no  output  links,  it  is  removed  from  this  list  and 
operation  of  the  algorithm  is  complete.  For  this  run  of  the  algorithm,  the 
earliest  arrival  times  of  the  message  at  N^,  N2,  and  N3  were  1,  2,  and  4, 
respectively. 

For  Monte  Carlo  run  2*  a coin  was  flipped  to  see  if  N2  was  destroyed  at 
time  3,  and  as  before  the  node  was  not  destroyed.  The  iterations  taken 
by  the  shortest  path  algorithm  are  given  in  Table  4.4 — iterations  0 and  1 were 
Identical  to  run  1.  During  iteration  2,  N2  was  considered  and  d23  was  selected  in 
accordance  with  the  distribution  given  in  Table  4.3.  A random  number  was 
generated  and  found  to  lie  below  0.667;  hence  d23  (2)  was  set  equal  to  1.  Since 

k*(2)  +d23-2+l<4-  k*(3) 
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TABLE  4.1 


Monte  Carlo  Run  1 of  the  Example 


ITERATION 

0 1 

■ 

2 

3 

CONSIDERA- 
TION LIST 

LEVEL 

1 

N1 

7 

1 

N2 

3 

■ 

■ 

B 

■ 

N3 

N3 

ARRIVAL 
LIST  (k*) 

NODE 

1 

1 

1 

1 

2 

00 

2 

2 

2 

3 

00 

B 

fl 

■ 
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TABLE  4.2 

Distributions  of  d^  and  d^ 

(First  Subscript  is  Sending  Node,  Second  Subscript  is  the  Link) 


J 

1 


3 


£ 

Prob(d^*£) 

E Pr 
m=l 

1 

0.667 

0.667 

2 

0.222 

0.889 

3 

0.074 

0.963 

4 

0.006 

0.969 

5 

0.005 

0.974 

• 

• 

• 

• 

• 

• 

• 

• 

• 

3 

0.667 

0.667 

6 

0.167 

0.834 

9 

0.083 

0.917 

• 

• 

• 

• 

• 

• 

• 

• 

• 

TABLE  4. 

3 

Distributions 

of  d23 
£ 

£ 

Prob (d23«L) 

E Pi 

m“l 

1 

0.667 

0.667 

2 

0.222 

0.889 

3 

0.018 

0.907 

4 

0.016 

0.923 

5 

0.0613 

0.939 

• 

• 

• 

• 

• 

• 

U 
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TABLE  4.4 


Monte  Carlo  Run  2 of  the  Example 

i 


< 


was  moved  from  level  4 to  level  3 on  the  consideration  list  and  k*(3) 
was  set  equal  to  3.  At  iteration  3 only  was  on  the  consideration  list. 
Since  it  has  no  links  leaving  it , it  was  deleted  and  this  run  was  over — 
the  earliest  times  of  arrival  of  the  message  at  N2  and  N3  were  1,  2, 
and  3,  respectively. 

In  this  simple  example  no  more  than  one  node  appeared  in  each  level 
of  the  consideration  list.  For  more  complicated  examples  each  level  may 
contain  many  nodes  at  any  given  iteration  of  the  algorithm.  The  order  of 
nodes  in  a given  level  is  immaterial,  but  placing  the  node  at  the  proper 
level  is  required  for  proper  operation  of  the  algorithm. 
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V. 


DETERMINATION  OF  STRATEGIC  C 


2 


SYSTEM  EFFECTIVENESS 


The  previous  three  sections  discussed  Monte  Carlo  simulation  of  the 
2 

Strategic  C Network  mission  of  delivering  attack  orders.  The  present  section 
addresses  three  questions:  What  are  useful  measures  of  system  effectiveness? 

How  many  Monte  Carlo  runs  are  required  to  accurately  determine  these  measures 
of  system  effectiveness?  How  can  the  number  of.  Monte  Carlo  runs  required  be 
reduced? 

Measures  of  System  Effectiveness 

It  is  convenient  to  divide  the  special  nodes  of  the  network  into  a 
number  of  classes;  for  example,  all  Mlnuteman  nodes  might  be  one  class,  all 
Submarine  nodes  a second  class,  and  all  Bomber  nodes  a third  class.  It  is  also 
appropriate  to  assign  to  each  node  in  a class  a time  varying  value  for  receipt 
of  the  EAM;  for  example,  the  Mlnuteman  nodes  might  suffer  attrition  with 
time  so  that  the  value  of  receiving  the  EAM  diminishes  with  time — in  this  case 
a suitable  value  for  receipt  of  the  EAM  is  the  number  of  Mlnuteman  Missiles 
remaining  at  the  time  of  message  receipt.  Given  the  classes  and  their  values 
for  receipt  of  the  EAM  at  various  times,  the  total  value  of  EAM  receipt  derived 
by  members  of  a given  class  at  a given  time  is  readily  determined  for  each 
Monte  Carlo  run.  With  these  results  for  a number  of  Monte  Carlo  runs,  several 
interesting  measures  of  system  performance  may  be  determined;  in  particular, 

time  period  histograms  of  value  and  arrival  time  distributions  are  of 
interest. 

Time  Period  Histograms.  For  this  measure  a particular  time  and  set  of  classes 
are  selected  and  a histogram  of  the  total  value  of  EAM  receipt  by  nodes  by  the 
specified  time  is  generated  (i.e.,  probability  of  obtaining  various  values 
where  value  is  appropriately  quantized).  For  example,  suppose  the  class  of  all 
Mlnuteman  nodes  is  selected  and  the  value  of  receiving  the  EAM  is  one  at  all 
nodes;  then  the  histogram  would  correspond  to  the  number  of  Mlnuteman  nodes 
receiving  the  EAM  by  a particular  time.  In  general,  several  families  of 
histograms — each  family  for  a given  set  of  classes,  each  member  of  the  family 
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for  a given  time — may  be  of  interest.  Some  examples  of  time  period  histograms 
are  presented  in  Section  VI. 

Arrival  Time  Distributions.  Consider  a set  of  classes  C, , C0,  ....  C^;  assume 
a set  of  minimum  values  V^j^;  Vj^,  ...»  V^^;  and  let  v^(t),  v2(t),  ...»  vn(t) 
be  the  values  for  EAM  receipt  by  time  t for  each  of  the  classes.  The  arrival 
time  distribution  is  the  cumulative  distribution  of  the  probability  that 
vl<T)  > VMjN,  v2(t)  > VMjN,  ...»  vn  (T)  > VM^  at  or  before  time  T For  example 
let  be  the  Minuteman  nodes,  C2  be  the  Submarine  nodes,  and  be  the  bomber 
nodes  and  suppose  it  is  desired  to  know  how  soon  n^  Minuteman  nodes,  n2 

Submarine  nodes,  and  n bomber  nodes  receive  the  message.  In  this  case,  the  value 

5 1 
of  receiving  the  EAM  at  each  special  node  is  set  at  unity  and  VMIN  * n^, 

V 2 - n_,  and  - n_.  The  result  will  be  the  cumulative  distribution  of 

MIN  2 MIN  3 

when  n^  Minuteman  nodes,  n2  Submarine  nodes,  and  n^  bomber  nodes  receive  the 
message.  Some  example  arrival  time  distributions  are  presented  in  Section  VI. 

Performance  Indices.  A single  measure  of  network  performance  is  the  expected 
total  value  of  special  nodes  which  receive  the  EAM  prior  to  specified  times 
of  interest.  For  example,  suppose  the  classes  are  the  three  TRIAD  forces 
as  above  and  the  value  of  EAM  receipt  is  proportional  to  the  size  of  the 
surviving  force  at  the  time  of  message  receipt  as  weighted  by  the  relative 
value  of  each  type  of  force.  The  resulting  performance  index  is  the  expected 
total  value  of  the  forces  receiving  the  EAM.  Again  an  example  of  a performance 
index  is  given  in  Section  VI. 

Statistical  Analysis  By  Monte  Carlo  Simulation 

It  is  clear  that  in  attempting  to  compute  a performance  index  one  is 
attempting  to  compute  the  expected  value  of  a random  variable.  The  same  is 
true  when  one  is  attempting  to  compute  a given  point  on  a time  period  histo- 
gram or  an  arrival  time  distribution.  For  example,  consider  the  problem  of 
computing  the  probability  of  deriving  a given  quantized  EAM  receipt  value  V 
on  a time  period  histogram.  Consider  the  random  variable  that  has  value  one 
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if  the  quantized  EAM  receipt  value  is  V and  zero  otherwise.  Then  the  proba- 
bility of  V is  the  expectation  of  this  random  variable.  In  general,  two 
properties  of  a Monte  Carlo  simulation  are  of  interest  for  each  random  var- 
iable: the  sample  mean,  used  to  estimate  the  expectation  of  the  random 
variable,  and  the  sample  variance,  used  to  measure  the  accuracy  of  that 
estimate.  From  these  two  properties  measures  of  confidence  may  be  generated. 


The  Sample  Mean  and  Sample  Variance.  Let  X be  a random  variable  and  ^ the 
value  of  X for  the  ith  Monte  Carlo  run.  Then  the  sample  mean  X and  the  sample 
variance  a2  for  m Monte  Carlo  runs  are  given  by: 


X 


(5-1) 


o2  = — 7 E (X  - X)2  (5-2) 

“-1  i-1  1 

Suppose  X has  mean  p and  variance  E.  The  central  limit  theorem  implies 
that  the  distribution  of  X approaches  normal  with  mean  \i  and  variance  E/m  as 
m increases;  furthermore,  if  X is  distributed  normally, 


is  distributed  via  Student's  t with  m degrees  of  freedom.  But  for  sufficiently 
large  m (say  m > 100)  Student's  t can  be  replaced  by  the  unit  normal  dis- 
tribution (i.e.,  normal  with  mean  zero,  variance  one).  In  the  sequel  it 
is  assumed  that  m is  sufficiently  large  such  that  t as  defined  in  (5-3)  is 
distributed  with  a unit  normal  distribution. 


Confidence  Measures.  From  the  above  it  is  clear  that  means  of  estimating  u 
for  a given  random  variable  X are  desired  and  that  X provides  such  an  estimate. 
Naturally,  the  question  of  how  good  the  estimate  of  p is  arises.  Suppose  t* 
is  selected  such  that  Prob  (t  > t*)  - 0.025,  where  t has  a unit  normal  distribu- 
tion. Suppose  that  100  Monte  Carlo  simulations  are  performed.  Each  simulation 
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may  measure  a different  random  variable;  the  only  requirement  is  that  a sufficient 
number  of  runs  are  made  for  each  simulation  to  insure  that  t in  (5-3)  is  unit 
normal.  It  may  be  expected  that  in  95  of  these  simulations  that  -t*  < t < t*  or 


X - t*  /a  / m < p < X + t* 


/a/m  . 


(5—4) 


This  approach  to  confidence  was  originated  by  E.  S.  Pearson  ; the  interval 

[X  - t*  JjTm  , X + t*  J?  /m  ] is  known  as  the  95%  two-sided  confidence 

interval  (other  confidence  intervals  are  defined  in  a similar  manner).  If 

every  time  a Monte  Carlo  simulation  is  made,  a 95%  confidence  interval  is  used, 

then  95%  of  the  time  X will  be  in  the  95%  interval.  However,  note  that  for  a 

given  Monte  Carlo  run  the  probability  that  X lies  in  the  interval  is  not  0.95— 

3 

in  fact  this  probability  is  not  defined. 


An  alternate  approach  to  confidence  was  proposed  by  R.  A.  Fisher.  His 
approach  is  as  follows.  Assume  that  X and  a are  fixed  at  the  values  obtained 
in  a Monte  Carlo  simulation  and  that  t is  unit  normal;  the  distribution  of 
t induces  a distribution  on  p known  as  the  fiducial  distribution.  Now  select 
and  u2  such  that  the  fiducial  probability  that  p < is  0.025  and  the 
fiducial  probability  that  p > P2  is  0.025.  The  interval  [p^  P2]  is  known 


as  the  95%  fiducial  interval;  clearly  in  the  present  case  p^  * X - t*  VO  /m 
— pi 

and  P2  " X + t*  /a  /m  and  the  confidence  and  fiducial  intervals  are 

3 

identical — this  is  not  always  the  case  however. 


Improving  Monte  Carlo  Efficiency 

Using  the  above  techniques  one  can  determine  when  enough  Monte  Carlo  runs 
have  been  made  to  estimate  any  desired  property  to  any  desired  accuracy 
as  specified  by  the  95%  (or  any  other  %)  confidence  intervals.  In  the  remainder 
of  this  section  two  methods  are  discussed  for  reducing  the  number  of  Monte 
Carlo  runs  required  to  achieve  a given  accuracy:  the  use  of  exact  comparisons 

and  the  use  of  biased  random  numbers.  The  reader  interested  in  pursuing  this 
matter  further  is  referred  to  reference  4. 


T 


Use  of  Exact  Comparisons.  Suppose  that  it  is  desired  to  compare  two  (or  more) 
alternative  Strategic  C2  networks  and  only  differences  in  performance  as 
opposed  to  absolute  performance  is  of  interest.  For  example,  a baseline 
network  may  exist  and  an  improved  network  may  be  proposed;  it  may  then  be 
desired  to  determine  if  the  improved  network  yields  a significant  performance 
improvement  compared  to  the  baseline.  Consider  a single  measure  of  performance. 
If  separate  Monte  Carlo  simulations  are  made  for  each  network,  the  observed 
differences  in  performance  between  the  two  networks  are  due  to  two  factors: 
the  difference  between  the  random  numbers  used  in  similar  situations  in  each 
simulation,  and  the  inherent  difference  in  the  capability  of  the  two  networks. 
Only  the  last  factor  is  of  interest.  Differences  in  performance  due  to  the 
first  factor  may  be  lessened  as  follows : let  the  same  random  numbers  be 

chosen  for  the  same  event  for  the  same  Monte  Carlo  run  in  the  simulation 
of  each  network,  compare  the  two  networks  run  by  run,  and  find  the  mean  and 
variance  of  the  differences  in  performance.  When  two  separate  Monte  Carlo 
runs  are  made  for  the  two  networks  the  variance  of  the  difference  in  perfor- 
mance is  the  sum  of  the  variances  in  performance  for  the  two  simulations. 

This  variance  is  larger — and  in  general  much  larger — than  the  variance  of 
the  difference  in  performance  obtained  on  a run-by-run  comparison  with  identical 
random  numbers;  hence  the  number  of  Monte  Carlo  runs  required  to  detect 
significant  differences  in  performance  between  the  two  networks  will  in 
general  be  much  fewer  if  the  exact  comparison  technique  is  used.  To  perform 
an  exact  comparison  using  the  Dynamic  Network  Simulator  the  random  numbers 
used  in  each  Monte  Carlo  run  for  node  destruction,  node  delays,  and  link 
delays  must  be  the  same  for  corresponding  nodes  and  links  in  the  two  networks 
(wherever  such  correspondence  exists)  . 

Use  of  Biased  Random  Numbers.  Suppose  that  it  is  desired  to  see  how  a network 
will  perform  under  stressing  situations.  Such  situations  may  be  emphasized 
in  a Monte  Carlo  simulation  by  biasing  the  random  numbers.  Typically,  when 
this  is  done  more  variations  in  behavior  of  the  network  (or  other  system 
elements  being  simulated)  are  obtained  with  the  same  number  of  Monte  Carlo 
runs.  As  a result,  a better  distribution  of  outcomes  is  obtained  with  the  same 
number  of  Monte  Carlo  runs  (i.e.,  the  distribution  of  t is  more  nearly  normal). 
Biasing  the  random  numbers  in  the  Dynamic  Network  Simulator  may  be  simply 
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accomplished  by  picking  random  numbers  with  non-uniform  distributions  on  the 
unit  interval.  (Longer  node  and  link  delays  are  obtained  when  the  random 
numbers  are  larger;  this  may  be  done  by  squaring — or  cubing,  etc. — the  random 
units  and  subtracting  from  1,  thus  biasing  the  random  number  toward  1. 

Earlier  node  failures  are  obtained  when  the  random  numbers  are  smaller;  this 
may  be  done  by  squaring  the  random  units  thus  biasing  the  random  number  toward 
0.) 


When  the  random  numbers  are  biased  the  computation  of  sample  mean  and 
sample  variance  must  be  modified  to  account  for  the  fact  that  stressing 
situations  are  occurring  more  often  than  expected.  The  necessary  modification 
is  to  use  weighted  averages  in  computing  the  sample  mean  and  variance  where 
the  weighting  for  each  Monte  Carlo  run  is  the  ratio  of  its  probability  of 
occurrence  in  nature  to  the  probability  of  its  occurrence  in  the  simulation. 

In  the  Dynamic  Network  Simulator  this  weighting  factor  is  found  as  follows: 
let  r, ....  r_  be  all  the  random  numbers  required  by  the  simulation,  let 


m 


th 


Pj (r)  be  the  probability  density  for  the  i 
and  let  be  the  value  of  r^  used  on  the  i 


random  number  in  the  simulation, 


th 


run.  Assuming  r^  is  a uniformly 


distributed  random  variable,  the  weighting  factor  for  the  i 


tb 


m 

IT 

j-l 


1/pJ(Cij) 


run  is 

(5-5) 


The  weighted  sample  mean  and  variance  are: 


X 


1 m 
- I 

m i-1 


Xi  pi 


(5-6) 


a 


2 


1 

m-1 


m 


(5-7) 


Note  that  if  unbiased  random  numbers  are  used  ■ 1 and  (5-6)  and  (5-7) 
reduce  to  (5-1)  and  (5-2). 
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VI . AN  EXAMPLE 


clearer  and  to  illustrate  the  type  of  outputs 
To  make  the  previous  discussion  clearer,  ana 

. His  an  example  is  presented  in  this  section.  It  should  be 

Ltm  — u • -«  - — - 

he  drawn  from  the  results. 


jut£Uts  Of  the  M»«-wnrk  Status  Model 


Consider  a hypothetical  network  of  33  nodes  and  97  links.  On  the  basis  of 
..  . od_  burst  and  jammer  information,  the  Network  Status  Model 
appropriate  link,  de  hn « - J ^ ^ ^ through  _ly  node  status 

generates  outputs  th  » « dMtructlon  versus  tine,  and  link  status 

in  terms  of  cumulative  probabil  y satellite, 

;r  ss 

M .......  « »•  " 

, h 1 1nk  If  one  inspects  Tables  6-4  through  6-6  caretuliy 

receiving  node  for  each  link.  Hp«.  2 3 13  and  26  through 

one  will  note  that  no  links  either  leave  or  enter  nodes  2,  3,  . 

32 i these  nodes  are  assumed  to  be  inoperable  for  this  example. 

ThP.  outputs  of  the  Dynam^  Nprwork  Simulator 

The  Dynamic  Network  Simulator  re,uires  for  inputs  the  outputs  of  the  Network 
star!  Model  such  as  those  given  in  Tables  6-1  to  6-6.  Assume  for  our  hypo  he 
tics!  network  that  the  outputs  shown  in  Tables  6-1  to  6-6  have  been  generated  y 
the  Network  Status  Model  and  supplied  to  the  Dynamic  Network  Simulator.  Note 
tie  Ilir  nodes,  as  described  in  Table  6-3,  are  divided  into  classes  »d 
the  values  given  in  that  table  for  receiving  the  message  at  a given  time. 

Dynamic  Network  Simulator  produces  a variety  of  outputs.  For  t e exanp 
Illustrative  outputs  are  shown  in  Figures  6-1  to  6-5  for  100  Monte  Carlo  . 
Piglre  6-1  to  L are  time  period  histograms,  taken  at  60  minutes  after  EAM 
Injection  into  the  network,  of  the  value  of  EAM  receipt  for  the  bomber  no  e . 
the  submarine  nodes,  the  Minuteman  nodes,  and  all  special  nodes.  Figure  6-5 

” “ I , „i  „„Mnt the  data  shown  in  Tables  6-1 

‘ “62  w^de^for  the^xsmple  and  not  produced  by  the  Network  Status 

Model. 
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TABLE  6.1  PROBABILITY  OF  NODE  DESTRUCTION  FOR  EXAMPLE 


NODE 

TIME 

PK 

TIME 

PK 

TIME 

PK 

TIME 

PK 

4 

.20 

10 

.40 

20 

.60 

5 

.20 

10 

.40 

20 

.60 

6 

.20 

10 

■ 

20 

.60 

7 

.20 

10 

20 

.60 

8 

B||l 

.10 

10 

'1 

20 

.30 

30 

.40 

9 

.10 

10 

.20 

20 

.30 

30 

.40 

10 

.10 

10 

.20 

20 

.30 

30 

.40 

11 

.20 

; 12 

.20 

20 

.10 

20 

o 

CM 

• 

21 

1 

.10 

20 

.20 

22 

1 

.50 

24 

1 

.20 

Notes: 

1.  PK  is  the  cumulative  probability  that  the  node  has 
been  destroyed  at  the  corresponding  time  or  sooner. 

2.  Nodes  not  appearing  in  table  are  assumed  to  be 
indestruc table . 
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TABLE  6.3  HIGH  FREQUENCY  LINKS  IN  EXAMPLE 


LINK 

N1 

N2 

CLASS 

(TIME,  PROBABILITY  OF  LINK  AVAILABILITY) 

1 

17 

20 

3 

(0, . 9)  (10,. 3)  (20, .7)  (30, .4) 

2 

17 

21 

3 

(0, .7) 

3 

1 

22 

2 

(0, .6)  (10,. 9)  (20,. 1)  (30,. 5) 

4 

22 

24 

2 

(0, .7) 

5 

22 

7 

2 

(0, .7) 

6 

22 

4 

2 

(0, . 7) 

7 

22 

6 

2 

(0,  .7) 

8 

22 

5 

2 

(0, .7) 

9 

1 

1? 

3 

(0, .7) 

10 

12 

11 

3 

(0, .7) 

11 

3 

(0, . 7) 

12 

3 

(0, .7) 

13 

3 

(0, .7) 

14 

3 

(0, .7) 

15 

3 

(0, .7) 

16 

3 

(0, .7) 

17 

3 

(0, .7) 

18 

17 

3 

(0, .7) 

19 

17 

3 

(0, .7) 

20 

3 

(0, .7) 

21 

5 

3 

(0, .7) 

22 

5 

6 

3 

(0, .7) 

23 

6 

5 

3 

(0, .7) 

24 

5 

8 

3 

(0, .7) 

25 

8 

5 

3 

(0, .7) 

26 

5 

18 

3 

(0, .6)  (10,. 9)  (20,. 1)  (30,. 5) 

27 

8 

6 

3 

(0, .7) 

28 

6 

8 

3 

(0, .7) 

29 

8 

15 

3 

(0, .7) 

30 

8 

9 

3 

(0,  .7) 

31 

9 

8 

3 

(0, .7) 

32 

9 

14 

3 

(0, .9)  (10, .3)  (20,. 7)  (30, .4) 

33 

8 

10 

3 

(0, .7) 

34 

8 

3 

(0, .7) 

35 

6 

3 

(0, .7) 

36 

ii  6 

E J 

3 

(0, .7) 

37 

BU 

3 

(0,  .7) 

38 

u 

3 

(0, .7) 

Notes: 

1.  N1  is  sending  node;  N2  is  receiving  node. 

2.  Class  2 links  are  satellite  links. 

3.  Link  probability  of  availability  is  assumed  to  hold  from  the  corresponding 
time  to  the  next  larger  time. 

I 
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TABLE  6.4  LOW  FREQUENCY  LINKS  IN  EXAMPLE 


N1 

N2  1 

MESS 

(TIME.  PROBABILITY  OF  CHARACTER  ERROR) 

1 

16 

13 

(0, .050) 
(25,. 001) 

(5, .040) 

(10,. 030) 

(15,. 020) 

(20,. 010) 

1 

19 

13 

(0, .020) 
(25,. 020) 

(5,. 040) 
(30, .010) 

(10,. 050) 
(35,. 005) 

Q5,  .040) 
'v40,  .005) 

(20,. 030) 

1 

6 

13 

(0..001) 
(25,. 020) 

(5,. 020) 
(3C..001) 

(10,. 030) 

(15,. 030) 

(20,. 010) 

1 

7 

13 

(0, .050) 
(25,. 001) 

(5,. 040) 

(10,. 030) 

(15,. 020) 

(20,. 010) 

1 

5 

13 

(0, .020) 
(25,. 020) 

(5,. 040) 
(30,. 010) 

(10,. 050) 
(35,. 005) 

(15,. 040) 
(40,. 005) 

(20,. 030) 

1 

4 

13 

(0..001) 
(25,. 020) 

(5,. 020) 
(30,. 001) 

(10,. 030) 

(15,. 030) 

(20,. 010) 

1 

17 

13 

(0, .050) 
(25, .001) 

(5,. 040) 

(10, .030) 

(15,. 020) 

(20,. 010) 

1 

18 

13 

(0, .020) 
(25,. 020) 

(5, .040) 
(30,  .010) 

(10, .030) 
(35,. 005) 

(15,. 040) 
(40, .005) 

(20,. 030) 

1 

15 

13 

(0..001) 
(25,. 020) 

(5, .020) 
(30, .001) 

(10, .030) 

(15, .030) 

(20, .010) 

7 

16 

13 

(0, .050) 
(25,. 001) 

(5,. 040) 

(10,. 030) 

(15,. 020) 

(20,. 010) 

7 

19 

13 

(0, .020) 
(25,. 020) 

(5, .040) 
(30,. 010) 

(10,. 050) 
(35,. 005) 

(15,. 040) 
(40,. 005) 

(20,. 030) 

7 

14 

13 

(0..001) 
(25,. 020) 

(5,. 020) 
(30,. 001) 

(10,. 030) 

(15,. 030) 

(20,. 010) 

7 

15 

13 

(0,-050) 

(25,-001) 

(5,. 040) 

(10,. 030) 

(15,-020) 

(20, .010) 

7 

18 

13 

(0, .020) 
(25,. 020) 

(5,. 040) 
(30,. 010) 

(10,. 050) 
(35,. 005) 

(15,. 040) 
(40,. 005) 

(20,. 030) 

7 

17 

13 

(0..001) 
(25,. 020) 

(5,. 020) 
(30,. 001) 

(10,. 030) 

(15,. 030) 

(20, .010) 

7 

6 

13 

(0, .050) 
(25,. 001) 

(5,. 040) 

(10,. 030) 

(15,-020) 

(20,. 010) 

7 

5 

13 

(0, .020) 
(25,. 020) 

(5,. 040) 
(30, .010) 

(10,. 050) 
(35,. 005) 

(15,. 040) 
(40,. 005) 

(20,. 030) 

7 

4 

13 

(0..001) 
(25,. 020) 

(5, .020) 
(30,. 001) 

(10,. 030) 

(15,-030) 

(20,-010) 

6 

7 

13 

(0, .050) 
(25,. 001) 

(5, .040) 

(10,. 030) 

(15,. 020) 

(20,. 010) 

6 

16 

13 

(0, .020) 
(25,. 020) 

(5,. 040) 
(30, .010) 

(10,. 050) 
(35,. 005) 

(15,. 040) 
(40,. 005) 

(20,. 030) 

6 

19 

13 

(0..001) 
(25,. 020) 

(5,. 020) 
(30,. 001) 

(10,. 030) 

(15,. 030) 

(20,. 010) 

6 

14 

13 

(0, .050) 
(25, .001) 

(5,. 040) 

(10,. 030) 

(15,. 020) 

(20,. 010) 

6 

15 

13 

(0, .020) 
(25,. 040) 

(5,. 040) 
(30,. 020) 

(10,. 050) 
(35,. 005) 

(15, .040) 
(40,. 005) 

(20,. 030) 

! 6 

18 

13 

(0..001) 
(25,. 020) 

(5,. 020) 
(30,. 001) 

(10,. 030) 

(15,. 030) 

(20,. 010) 

j (0, .050)  (5,. 040)  (10,. 030)  (15,. 020)  (20,. 010) 


(25,. 001) 


6 


5 
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TABLE  6.4  LOW  FREQUENCY  LINKS  IN  EXAMPLE  (CONT.) 


LINK 

N1 

N2 

MESS 

(TIME,  PROBABILITY 

OF  CHARACTER  ERROR) 

64 

6 

17 

13 

(0, .020) 
(25,. 020) 

(5,-040) 
(30,. 010) 

(10,-050) 

(35,-005) 

(15,. 040) 
(40,. 005) 

(20,-030) 

65 

6 

4 

13 

(0,.001) 
(25,. 020) 

(5,-020) 
(30,. 001) 

(10,. 030) 

(15,. 030) 

(20, .010) 

66 

4 

7 

13 

(0, .050) 
(25,. 001) 

(5, .040) 

(10,. 030) 

(15,. 020) 

(20, .010) 

67 

4 

6 

13 

(0, .020) 
(25,. 020) 

(5,-040) 
(30, .010) 

(10,. 050) 
(35,. 005) 

(15,. 040) 
(40, .005) 

(20, .010) 

68 

4 

19 

13 

(0..001) 
(25,. 020) 

(5,. 020) 
(30,. 001) 

(10,. 030) 

(15, .030) 

(20, .010) 

69 

4 

16 

13 

(0, .050) 
(25, .001) 

(5,. 040) 

(10,. 030) 

(15,-020) 

(20,-010) 

70 

4 

14 

13 

(0, .020) 
(25,. 020) 

(5,  .040) 
(30,  .010) 

(10,. 050) 
(35,. 005) 

(15, .040) 
(40,-005) 

(20,. 030) 

71 

4 

15 

13 

(0..001) 
(25,. 020) 

(5, .020) 
(30, .001) 

(10,. 030) 

(15,-030) 

(20,. 030) 

72 

4 

18 

13 

(0, .050) 
(25,. 001) 

(5,. 040) 

(10,. 030) 

(15,-020) 

(20,. 010) 

73 

4 

17 

13 

(0, .020) 
(25,. 020) 

(5,. 040) 
(30,. 010) 

(10,. 050) 
(35,. 005) 

(15,-040) 

(40,-005) 

(20,. 020) 

74 

5 

6 

13 

(0..001) 
(25,. 020) 

(5,. 020) 
(30, .001) 

(10,. 030) 

(15,-030) 

(20,. 010) 

75 

5 

7 

13 

(0,-050) 
(25,. 001) 

(5,. 040) 

(10,. 030) 

(15,-020) 

(20,-010) 

76 

5 

4 

13 

(0, .020) 
(25,-020) 

(5,. 040) 
(30,. 010) 

(10,. 050) 
(35,. 005) 

(15, .040) 
(40, .005) 

(20,-020) 

77 

5 

17 

13 

(0,-001) 

(25,-020) 

(5, .020) 
(30,. 001) 

(10,. 030) 

(15,. 060) 

(20,-010) 

78 

5 

18 

i 13 

(0,-050) 

(25,-001) 

(5, .040) 

(10,. 030) 

(15,-020) 

(20,. 010) 

79 

5 

15 

13 

(0,-020) 

(25,-020) 

(5,. 040) 
(30,. 010) 

(10,. 050) 
(35,. 005) 

(15, .040) 
(40,. 005) 

(20,. 020) 

80 

5 

14 

13 

(0,-050) 
(25, .001) 

(5,. 040) 

(10,. 030) 

(15,-020) 

(20,. 010) 

81 

5 

16 

13 

(0..001) 

(25,-020) 

(5,. 020) 
(30,-001) 

(10,. 030) 

(15,-030) 

(20,-010) 

82 

5 

19 

13 

(0,-050) 

(25,-001) 

(5,-040) 

(10,. 030) 

(15,-020) 

(20, .010) 

83 

4 

5 

13 

(0,-040) 

(25,-001) 

(5,-040) 

(10,. 030) 

(15,-020) 

(20,-010) 

84 

24 

33 

13 

(0,-020) 

(25,-020) 

(5,-030) 

(30,-010) 

(10,. 050) 
(35,. 005) 

(15,-040) 

(40,-005) 

(20,-020) 

Notes: 

1.  N1  is  sending  node;  N2  is  receiving  node. 

2.  All  low  frequency  llnkB  are  in  class  4. 

3.  Link  probability  of  character  error  is  assumed  to  hold 
from  corresponding  time  to  next  larger  time. 

4.  "MESS"  is  the  time  required  for  a single  transmittal  of  the  message. 
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(TIME,  PROBABILITY  OF  AVAILABILITY) 


(0, .7) 

(0..7) 

(0,  .9)  (10, .3)  (2C , .7 ) (30, .4) 

(0, .7) 

(0, .7) 

(0, .7) 

(0, .7) 

(0,  .7) 

(0, .7) 

(0, . 7) 

(0, .7) 

(0, .7) 

(0, .6)  (10, .9)  (20, .1)  (30, .5) 


1.  N1  is  sending  node;  N2  is  receiving  node. 

2.  All  special  links  are  class  1 links. 

3.  Time  Delay  is  the  time  delay  after  EAM  arrival  before  the  link  becomes 
available . 

4.  Link  probability  is  assumed  to  hold  from  the  corresponding  time  to  the 
next  larger  time. 


* 
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TABLE  6.6  SPECIAL  NODES  AND  THEIR  PROPERTIES 


IODE 

CLS 

(TIME.  VALUE) 

PAIRS 

23 

1 

(0,499.00) 

(10,450.00) 

(20,375.00) 

(30,250.00) 

25 

1 

(0,499.00) 

(10,450.00) 

(20,375.00) 

(30,250.00) 

33 

2 

(0,999.00) 

14 

3 

(0,200.00) 

(10,100.00) 

(20,40.00) 

( 30,30.00) 

15 

3 

(0,150.00) 

(10,70.00) 

(15,40.00) 

(20,30.00) 

16 

3 

(0,200.00) 

(10,100.00) 

(20,40.00) 

(30,30.00) 

17 

3 

(0,150.00) 

(10,70.00) 

(15,40.00) 

(20,30.00) 

18 

3 

(0,200.00) 

(10,100.00) 

(20,40.00) 

(30,30.00) 

19 

3 

(0,150.00) 

(10,70.00) 

(15,40.00) 

(20,30.00) 

Notes: 

1.  Class  1 nodes  are  bomber  nodes;  class  2 nodes  are 
submarine  nodes;  class  3 nodes  are  Minuteman  nodes. 

2.  Value  of  the  EAM  to  a node  Is  assumed  constant  from 
the  associated  time  to  the  next  larger  time. 

t 
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Probability 


100  200  300  4 DO  500  600  700 

Units  of  Value  for  Bomber  Nodes 


FIOURF.  6.1  TIME  PERIOD  HISTOGRAM  OF  VALUE  ACCRUED  BY  BOMBER  NODES  FOR 
HAVING  RECEIVED  THE  EAM  AT  OR  PRIOR  TO  t-60  MINUTES 


Probability 


300  400  500  600  700  800  900  1000 

Units  of  Value  for  Submarine  Nodes 


FIGURE  6.2  TIME  PERIOD  HISTOGRAM  OF  VALUE  ACCRUED  BY  SUBMARINE  NODES 
FOP  HAVIN'  RECEIVED  THE  EAM  AT  OR  PRIOR  TO  r-60  MINUTES 


Probability 


Probability 


< 

i 


100  200  300  400  500  600 

Units  of  Value  for  Minutemen  Nodes 


FIGURE  6.3  TIME  PERIOD  HISTOGRAM  OF  VALUE  ACCRUED  BY  MINUTEMAN  NODES 
FOR  HAVING  RECEIVED  THE  EAM  AT  OR  PRIOR  TO  t-60  MINUTES 


.0  _ 

.8  - 


250  500  750  1000  1250  1500  1750  2000  2250 


Units  of  Value  for  all  Tbree  Classes  Combined 


FIGURE  6.4  TIME  PERIOD  HISTOGRAM  OF  VALUE  ACCRUED  BY  ALL  SPECIAL  NODES 
FOR  HAVING  RECEIVED  THE  EAM  AT  OR  PRIOR  TO  t-60  MINUTES 
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plots  the  cumulative  probability  that  certain  combinations  of  value  are  exceeded 
versus  time. 

Suppose  that  the  final  time  of  interest  is  60  minutes  (i.e.,  value  of 
receiving  the  message  after  that  time  is  zero).  Then  Figure  6-4  is  a histogram 
of  the  performance  index  (i.e.,  the  total  value  of  the  nodes  that  receive  the 
EAM).  The  expected  performance  is  J = 860. 


VII.  EXTENSIONS  OF  THE  PRESENT  MODEL 


This  report  has  described  two  models  --  the  Network  Status  Model  and  the 

Dynamic  Network  Simulator  that  can  be  used  to  simulate  in  an  efficient  manner 

o 

the  behavior  of  a Strategic  C network  delivering  attack  orders.  The  Network 

2 

Status  Model  should  be  applicable  with  little  change  to  other  Strategic  C 
problems  such  as  the  Attack  Warning  and  Battle  Management  problems.  On  the 
other  hand,  the  Dynamic  Network  Simulator  in  its  present  form  is  applicable 

only  to  deliverv  of  attack  orders.  The  heart  of  the  Dynamic  Network  Simulator 
is  a version  of  the  shortest  path  algorithm.  Direct  use  of  the  SPA  is  possible 

because  it  has  been  assumed  that  onlv  one  message  is  transmitted  by  the  network 

and  that  message  is  transmitted  to  all  nodes.  In  the  Attack  Warning  and  Battle 

Management  problems  more  than  one  message  is  sent  over  the  network  and  each 

message  is  sent  to  only  one  of  a few  nodes.  For  these  problems  routing  and 

priority  of  messages  become  important  and  the  DNS  must  be  modified  to  simulate 

2 

the  behavior  of  the  Strategic  C Network  employed. 

The  Shortest  Path  Algorithm  Viewed  as  a Discrete  Simulation 

In  order  to  extend  the  Dynamic  Network  Simulator  to  address  the  Attack 
Warning  and  Battle  Management  problems  it  is  helpful  to  view  it  from  a different 
point  of  view.  A discrete  event  simulation  is  one  in  which  events  are  time 
ordered  and  performed  in  time  sequence;  any  event  may  create  new  future  events, 
cancel  future  events,  or  change  the  time  of  future  events.  The  basis  of 
a discrete  event  simulation  is  the  event  list  containing  the  events  to  be 
performed;  the  simulation  proceeds  by  performing  the  earliest  event  on  the 
list.  Now,  if  events  are  defined  to  be  arrival  of  the  message  at  a node 
then  it  is  clear  that  the  SPA  as  used  in  the  Dynamic  Network  Simulator  is  a 
discrete  event  simulation  of  a Strategic  C Network  and  that  the  consideration 
list  is  the  event  list  for  the  simulation. 

Communication  networks  have  been  modelled  by  discrete  event  simulations 
in  the  past  (see  for  example  Ref.  5);  however,  the  simulation  in  Ref.  5 required 
significantly  more  computer  time  than  that  of  the  simulation  described  in 
this  report.5  The  reasons  for  the  relative  computer  time  requirements  of 
the  two  simulations  are  threefold;  computers  are  now  faster,  the  previous 
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simulation  used  more  than  the  minimum  number  of  events  that  are  necessary,  and 
the  previous  simulation  did  not  use  the  efficient  organization  provided  by  the 
SPA.  The  latter  two  reasons  relate  to  the  efficiency  of  processing  the  event 
list,  which  in  turn  depends  upon  how  the  event  list  is  ordered.  If  the  event  list 
contains  events  listed  in  order  of  their  creation  then  it  must  be  searched  each 
time  an  event  is  created,  deleted  or  modified.  Such  searches  will  typically  be 
linear  — hence  the  search  time  is  proportional  to  the  length  of  the  list.  How- 
ever, if  a binary  list  is  used  a binary  search  (wnich  has  search  time  propor- 
tional to  the  logarithm  of  list  length)  is  possible  when  placing  new  events  in 
the  list  but  not  when  deleting  or  modifying  existing  events  on  the  list. 

I | 

A discrete  event  simulation  based  upon  the  SPA  as  described  in  this  report 
is  superior  to  standard  discrete  event  simulations  in  two  respects.  First, 
message  arrivals  are  the  only  events  considered;  any  event  — such  as  weapon 
launch  — that  is  dependent  on  message  arrival  can  readily  be  modelled  as  part 
of  the  message  arrival  event,  and  thus  the  number  of  events  that  must  be  con- 
sidered is  minimized.  Secondly,  by  quantizing  time  and  using  the  node  list, 
the  consideration  list  need  not  be  completely  searched.  To  see  this,  view  the 
consideration  list  as  a two  dimensional  array  where  one  dimension  is  quantized 
time  and  the  second  dimension  is  all  events  that  occur  at  that  time  in  any  order 
(e.g. , all  nodes  that  receive  a message  at  that  time  in  any  order).  It  is  clear 
that  if  the  information  contained  in  the  node  list  is  used,  creating,  modifying 
or  deleting  an  event  involves  only  a search  on  the  second  dimension  of  this 
two  dimensional  array.  Since  the  number  of  nodes  receiving  a message  at  a 
given  time  is  small,  the  required  search  time  is  small.  The  major  drawback  to 
this  approach  is  the  excessive  storage  that  would  be  required  if  very  small 
time  quantization  levels  were  required. 

Message  Routing  and  Priority 

The  discrete  event  simulation  developed  for  the  Delivery  of  Attack  Orders 
Problem  may  be  readily  modified  to  treat  the  Attack  Warning  and  Battle  Manage- 
ment Problems.  When  so  modified, its  heart  is  no  longer  a SPA  but  it  retains 
the  list  processing  efficiency  of  the  SPA.  The  major  required  modifications 
relate  to  incorporation  of  message  priority  and  routing  logic.  The  capacity  of 
links  is  limited  to  a small  number  of  messages  at  one  time  over  a high 
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frequency  link  and  one  message  at  a time  over  a low  frequency  link.  If  more 
messages  must  be  sent  over  a link  than  its  capacity  allows,  the  sequence  in 
which  messages  are  sent  must  be  determined.  For  a given  message, if  no  direct 
link  exists  between  the  node  under  consideration  and  the  ultimate  destination 
of  the  message*  then  the  next  node  to  which  the  message  is  sent  must  be 
selected.  Generation  of  the  logic  for  routing  messages  is  a task  that  is  made 
no  easier  by  the  fact  that  the  sender  may  be  unaware  of  the  status  of  the  links 
he  has  available  for  sending  messages.  An  important  use  of  a modified  Dynamic 
Network  Simulator  would  be  in  comparing  the  efficiency  of  various  existing  and 
proposed  routing  logics. 
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APPENDIX  A.  NODE  DESTRUCTION  ALGORITHM 


For  nodes  located  on  the  ground  and  for  ground  level  nuclear  bursts,  the 
node  kill  probabilities  are  computed  by  a standard  nuclear  weapons  effect  algorithm 
in  the  Network  Status  Model.  This  model  sets  the  probability  of  destruction  of 
the  node  equal  to  the  circular  coverage  function  P(R,r,a)  which  is  defined  to  be: 

P(R,r,a)  ■ J N(x,y)dx  dy, 
where  D 

N(x,y)  ■ 2-dimensional  normal  distribution  function 
for  mean  0 and  covariance  1, 

D - | (x,y) : (u-x)2  + (v-y)2  < (R/a)2,  u2+  v2  - (r/o)2  } . 

Note  that  from  the  axial  symmetry  of  N(x,y),  the  integral  is  independent  of  the 

2 2 2 

value  of  u and  v chosen  provided  u + v = (r/o)  . 


Very  briefly  the  model  is  as  follows.  The  parameter  r is  taken  to  be  the 
ground  distance  between  the  node  and  the  burst,  while  R is  computed  from  the  VN 
hardness  of  the  target  and  the  warhead  size.  The  parameter  a is  found  from 


where  aD  is  a known  function  of  R which  depends  upon  whether  the  target  is  over- 
pressure or  dynamic  pressure  sensitive  (indicated  in  the  VN  number  by  a P or  Q) , 

and  a is  given  by 
£ 


Og  * (CEP  of  warhead)  /1.1774 

It  is  assumed  that  nodes  are  point  targets;  if  this  is  not  true  a term  a2  must 

2 " 
be  added  to  the  equation  of  a in  order  to  reflect  the  size  of  the  node. 

The  computation  of  R proceeds  as  follows . From  the  VN  number  and  warhead 
size,  the  overpressure  or  dynamic  pressure  at  which  the  node  will  be  killed  50% 
of  the  time  is  computed  from  known  equations. ^ R is  the  ground  radius  from  the 
burst  point  at  which  this  overpressure  or  dynamic  pressure  occurs.  In  the  model 
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R is  found  from  an  approximate  equation  fitted  to  tabular  data.  Two  fits  are 
possible:  one  for  ground  level  bursts  and  one  for  bursts  at  optimum  height 
above  ground.  The  input  nuclear  weapon  laydown  tape  must  indicate  for  each 
ground  burst  whether  it  is  an  optimum  or  zero  height  burst. 


APPENDIX  B.  USE  OF  SHORTEST  PATH  ALGORITHM  TO 


DETERMINE  STATIC  CONNECTIVITY 


The  shortest  path  algorithm  presented  in  the  body  of  this  report  may  be 
used  to  predict  static  connectivity  of  source  to  sink  nodes  for  networks  with 
random  failures  of  nodes  and  links.  Let 


probability  node  i is  not  operable  at  some  specified  time 
probability  link  ij  is  not  operable  at  some  specified  time. 


Before  each  call  of  the  shortest  path  algorithm  during  each  Monte  Carlo  run,  a 
random  unit  r is  selected  for  each  node  i.  If  r < q then  the  time  of  node 
destruction  t°  is  set  to  0;  otherwise  t^  = 00 . Each  time  a link  is  considered 
in  the  shortest  path  algorithm,  a random  unit  r^  is  selected.  If  r^  < P^  then 
the  link  delay  d^  is  set  to  °°;  otherwise  d^  * 0.  Nodes  with  k*(i)  = 0 at  the 
end  of  each  Monte  Carlo  run  are  connected  to  the  source;  nodes  with  k*(i)  « 00 
are  not  connected  to  the  source. 


TVo  examples  were  run  to  compare  the  shortest  path  algorithm  with  a standard 
connectivity  algorithm  (SCA)1.  The  example  illustrated  in  Fig.  B.l  is  sufficiently 
simple  that  the  connectivities  with  the  source  node  may  be  easily  computed  by 
hand  — these  values  are  given  in  Table  B.l.  Two  different  results  are  given  for 
nodes  7 and  8 depending  upon  whether  links  7-8  and  8-9  are  directed  or  not.  If  p 
is  the  probability  of  being  disconnected, then  the  standard  deviation  in  a Monte 
Carlo  determination  of  p is  p(l-p)/N  where  N is  the  number  of  Monte  Carlo  runs. 
Computed  standard  deviations  for  N * 900  are  also  given  in  Table  B.l.  Monte  Carlo 
determinations  of  p for  N ■ 900  from  both  the  standard  connectivity  algorithm 
(SCA)  and  shortest  path  algorithm  (SPA)  are  given  in  Table  B.l,  as  well  as  the 
errors  in  these  estimates  normalized  by  the  appropriate  standard  deviations. 

Both  sets  of  answers  closely  approximate  the  true  probabilities  of  dis- 
connection. The  average  standard  deviations  are  biased  negatively,  but  this 
is  probably  caused  by  the  correlation  between  connectivities  as  the  following 
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LINK  FAILURE  PROBABILITIES 
-0.2  FOR  ALL  LINKS 


FIG.  B-l  AN  EXAMPLE  NETWORK 
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TABLE  B-l 


RESULTS  FOR  EXAMPLE  OF  FIG.  B-l 


R 

PROBABILITY  OF 
DISCONNECTION* 

STANDARD 

DEVIATION 

ESTIMATED  PROBABILITY  OF 
DISCONNECTION 
SCA  SPA 

ERROR /STANDARD  DEVIATIOl 
SCA  SPA 

2 

0.424 

0.016 

0.400 

0.403 

-1.5 

nm 

3 

0.539 

0.017 

0.490 

0.507 

-2.59 

4 

0.631 

0.016 

0.632 

0.624 

0.06 

5 

0.424 

0.016 

0.399 

0.422 

-1.57 

-0.12 

6 

0.539 

0.017 

0.525 

0.531 

-0.80 

-0.47 

7 

0.429 

0.016 

0.387 

-2.65 

7** 

0.488 

0.016 

0.474 

-0.87 

8 

0.406 

0.016 

0.377 

-1.86 

8** 

0.631 

0.016 

0.632 

+0.06 

9 

0.301 

0.015 

0.281 

0.270 

-1.33 

-2.07 

Ave.  -1.53  -0.87 


* Disconnection  here  implies  Node  1 and  Node  j . 
**  Link  Directed. 
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analysis  shows.  In  the  SPA  Monte  Carlo  runs  it  was  found  that  the  source  node 
failed  19.1%  of  the  time  as  compared  to  a computed  value  of  20.0%.  The  normalized 
error  is  -0.69  standard  deviations;  however,  since  failure  of  the  source  node 
implies  no  connectivity,  this  result  biases  all  the  connectivities  on  the  nega- 
tive side. 


A more  complicated  example  shown  in  Table  B.2  was  considered.  In  this  example 
there  are  34  nodes  (only  31  of  which  are  operable),  each  with  a probability  of  0.1 
of  failure,  and  91  links,  each  with  a probability  0.2  of  failure.  Both  the  SCA 
and  the  SPA  programs  were  used  to  determine  the  probabilities  that  nodes  were  con- 
nected with  the  source  node  by  making  100  Monte  Carlo  runs;  results  are  given  in 
Table  B.3,  in  which  X^,  X2»  and  X^  are  the  SCA  (undirected),  SPA  undirected,  and 
SPA  directed  probabilities  of  connection.  (With  the  SPA  program  connectivity  was 
computed  both  for  links  directed  as  in  Table  B.2  and  for  bi-directional  links  by 
assuming  a uni-directional  link  each  way  for  each  link  of  Table  B.2.)  Standard 
deviations  for  the  network  with  bi-directional  links  were  computed  by  the  SCA 
program.  Table  B.3  also  gives  the  error  between  the  bi-directional  results  from 
the  two  programs  normalized  by  J 2 times  the  standard  deviation  (because  each  is 
in  error  from  the  true  value  with  a covariance  equal  to  the  standard  deviation 
squared)  Again  the  average  normalized  difference  is  biased,  largely  because  in 
the  SCA  program  the  source  node  1 failed  only  9 times  and  in  the  SPA  program  it 
failed  only  4 times  instead  of  the  expected  10  times.  To  remove  this  effect 


0.90 

0.91 


V 


0.90  x 
0.96  2 


were  computed;  the  results  as  well  as  the  normalized  differences  between  X^  and 
X,.  are  displayed  in  Table  B.3.  The  average  normalized  difference  between  the  two 
for  bi-directional  links  is  now  0.07  standard  deviations,  indicating  that  two  pro- 
grams give  essentially  the  same  results.  The  time  taken  on  a UNIVAC  1108  computer 
to  obtain  the  results  of  Table  B.3  via  the  SPA  was  1.5  seconds  to  find  the  con- 
nectivities with  directed  links  and  2.8  seconds  to  find  the  connectivities  with 
undirected  links. 
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APPENDIX  C.  ORGANIZING  THE  DATA  BASE 


2 

In  typical  Strategic  C networks  of  interest  for  analysis  there  could  be  as 

many  as  700  links  and  300  nodes.  Each  of  these  elements  has  its  own  particular 

set  of  characteristics  as  well  as  a generic  set  of  characteristics  representative 

of  its  element  type.  This  represents  a sizable  data  base  that  must  be  properly 

organized  and  managed  as  a first  step  in  understanding  and  improving  a Strategic 

C system  as  a whole.  The  need  for  such  structuring  becomes  even  more  apparent 

when  trying  to  assess  the  effectiveness  of  a number  of  postulated  improvements 

2 

in  a network  and  its  components.  The  objective  of  organizing  the  Strategic  C data 
base  to  facilitate  effectiveness  studies  may  be  met  by  a computerized  data 
management  and  information  retrieval  system.  Development  of  this  system  would 
require  the  accomplishment  of  three  identifiable  sub-tasks:  designing  suitable 

formats  for  data  storage  on  various  media;  creating  the  software  for  data  input, 
retrieval,  and  editing;  and  obtaining  the  baseline  data  and  implementing  the 
system.  Thereafter  the  data  base  must  be  maintained  in  a useful  state  by  purging 
outdated  entries,  adding  new  entries,  and  modifying  spurious  entries. 


Data  Requirement 

Computer  processable  records  may  be  distinguished  on  the  basis  of  their 
length  attribute:  fixed  or  variable.  The  length,  is  measured  in  bytes  (which 

is  synonomouA with  characters  in  this  context),  and  'variable'  implies  that  each 


record  in  the  set  may  have  a different  length  while  'fixed'  implies  that  each 
record  has  the  ,^ame  length.  Though  variable  records  offer  certain  economies  of 
storage  they  arA clumsy  to  manage  in  data  retrieval  and  sorting  applications;  for 
this  reason  fixecn  records  are  preferred  in  the  present  application.  Each  fixed 
record  may  be  subdivided  into  fields  of  contiguous  strings  of  bytes  with  fixed 
starting  positions  \\o  record  individual  data  items.  The  length  of  these  fields 


may  vary 


within  a r^Aord,  but  is  the  same  for  all  records  within  the  data  set. 


Since  the  field  must  Contain  values  of  variables,  its  length  must  be  adequate  to 
encompass  the  entire  range  of  variation  of  the  variable  to  which  it  is  assigned. 
Thus  four  bytes  must  beyprovided  to  record  frequencies  in  the  range  0000  to 
1000  MHZ  with  a resolutidX  of  1 MHz. 


2 

The  most  convenient  organization  of  the  C data  base  is  to  generate  two  data 
sets  of  fixed  records;  one  data  set  being  allocated  to  link  data  and  the  other 
to  node  data.  A record  in  either  data  set  can  then  be  used  to  record  all  pertinent 
information  about  a particular  link  or  node.  Schematically  each  data  set  would 
then  have  the  following  organization: 


1st 

2nd 

3rd 

Data 

Data 

Data 

Item 

Item 

Item 

• • • • 

Element 

1 

1 1 

1 1 _ _ 

Element 

2 

! r~ 

I 1 

Element 

3 

i i 

1 1 ‘ 

• 

• 

• 

• 

• 

• 

• 

The  ordering  of  data  items  within  records  is  less  important  than  a determination  of 
what  those  items  should  be.  This  determination  must  be  made  on  the  basis  of  present 
and  fu'ure  availability  of  data  balanced  against  a judgement  of  what  data  is 

2 

desired;  thus  a two-pronged  effort  is  indicated.  Ideally,  elements  of  the  C 
system  must  be  analyzed  to  extract  the  relevant  characteristics  of  each  element 
that  affect  overall  system  performance,  regardless  of  whether  or  not  this  data  is 
available.  On  the  other  hand,  a literature  survey  must  be  conducted  to  ascertain 
the  actual  and  projected  availability  of  data. 

An  effort  such  as  that  outlined  above  can  isolate  the  generic  data  items,  but 
suffers  from  a lack  of  specificity  with  regard  to  particular  elements.  That  is, 
fields  within  the  record  nust  be  reserved  for  status  information  such  as:  a flag 

to  indicate  if  that  element  is  operative;  indicators  of  significant  departures 
from  operating  norms  for  that  type  of  element;  etc.  Equally  important  is  inter- 
connectivity information  — nodes  for  example  must  have  pointers  that  indicate 
other  nodes  that  address  it,  as  well  as  pointers  to  nodes  that  it  addresses.  The 
links  over  which  this  communication  occurs  must  also  be  recorded  as  pointers  from 
the  node  data  set  to  the  link  data  set  and  vice  versa. 
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Software  Requirements 


All  of  the  data  discussed  above  has  the  added  dimension  of  temporal  var- 

2 

iation:  nodes  and  links  are  added  or  deleted  from  the  C system,  link  charact- 

eristics may  suffer  diurnal  and  annual  variation,  nodes  are  upgraded  by  hard- 
ware modifications  and  downgraded  by  hardware  malfunctions.  These  variations 
imply  that  fast,  convenient  means  for  modifying  the  data  base  and  for  displaying 
status  information  must  be  sought.  This  in  turn  outlines  the  requirements  for 
several  desirable  computer  programs  which  arc  discussed  below. 

Data  Input  and  Validation  Program 

There  is  a significant  step  between  possessing  pertinent  data  and  transcribing 
it  on  its  ultimate  storage  medium  (tape,  disc,  or  magnetic  drum)  in  the  required 
format.  Normally,  printed  data  is  transcribed  on  punched  cards  that  are  read  by 
computer  and  transferred  to  the  storage  device.  This  error  prone  process  is  com- 
plicated by  the  necessity  of  using  other  than  the  desired  final  format  on  the 
punched  cards.  The  effect  of  keypunching  errors  may  be  minimized  by  designing 
a simplified  format  for  the  punched  cards,  and  using  a program  to  reformat  the 
data  and  validate  it  at  the  same  time.  For  example,  the  validation  program 
could  reject  an  entry  for  a node  having  a transmit  frequency  of  10  Hz  connected 
to  a tropospheric  scatter  link. 

Data  Base  Update  Program 

This  program  would  provide  three  facilities  and  operate  on  the  stored  data 
base:  ' 

1)  Modification  - a facility  to  correct  spurious  entries  in  the  data  base, 
or  to  alter  data  items  that  have  changed  since  the  last  updating  period. 

2)  Insertion  - a facility  for  adding  new  elements  or  records  to  the 
existing  data  sets. 

3)  Deletion  - a facility  for  deleting  elements  or  records  from  the  existing 
data  sets. 
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This  program  must  be  designed  with  a simplified  language  for  performing  its 
various  functions  so  that  a minimum  of  training  would  be  required  to  use  it. 

At  the  same  time,  it  must  have  built  in  safeguards  to  prevent  inadvertent  destruc- 
tion of  valuable  data.  At  the  very  least  this  implies  that  the  update  program 
mjst  provide  an  automatic  backup  capability  to  preserve  the  data  sets  in  unedited 
form. 

Data  Retrieval  and  Display  Programs 

This  may  be  one  or  a collection  of  programs  to  interface  between  the 
storage  media  and  the  data  user.  An  interactive  system,  perhaps  utilizing  a 
cathode  ray  tube,  would  be  highly  desirable.  The  system  would  provide  a sim- 
plified retrieval  language  enabling  the  user  to  specify  his  needs.  Computer 
prompting  could  be  incorporated  to  lead  the  user  through  a decision  tree  of 
possible  data  alternatives.  These  programs  must  have  the  ability  to  sort  the 
data  on  various  keys,  retrieve  specific  data  items,  and  generate  simple  statistics. 
Some  example  commands  that  this  system  could  support  are: 

DISPLAY  STATUS  OF  LINK  230 

LIST  INPUTS  TO  NODE  570 

GIVE  PERCENT  OF  OPERATIONAL  VLF  LINKS 

LIST  ALL  SEQUENCES  OF  LINKS  BETWEEN  WASHINGTON  D.C.  AND  OMAHA 

To  increase  the  versatility  of  these  programs  there  must  be  provision  for 
various  data  display  options.  The  cathode  ray  tube  is  suitable  for  displaying 
items  of  transitory  interest  such  as  the  man-computer  dialogue,  and  to  display 
graphical  output  such  as  network  segments.  A hard-copy  facility  must  also  be 
provided  to  permanently  record  the  generated  graphs,  if  desired,  and  to  record 
voluminous  tabular  displays.  Careful  thought  must  be  invested  in  designing 
these  displays  to  insure  that  meaningful  data  is  highlighted,  and  that  less 
meaningful  data  is  suppressed  unless  specifically  requested. 

The  ground  work  for  an  effective  and  versatile  data  management  and  retrieval 
2 

system  for  Strategic  C has  been  outlined  above.  Implementation  issues  revolving 
mainly  around  questions  of  host  hardware  can  not  be  addressed.  It  is  axiomatic 
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that  the  larger  and  faster  the  host  computer,  the  greater  is  the  generality,  flexi- 
bility and  responsiveness  that  can  be  achieved.  In  order  of  increasing  preference, 
the  host  computer  could  provide  tape,  disc,  or  drum  storage,  with  combinations 
of  these  three  preferable  to  any  one  singly.  For  interactive  retrieval,  various 
CRT  or  typewriter  terminals  are  desirable,  and  a line  printer  for  hard  copy 
output  is  mandatory.  The  core  requirement  is  a strong  function  of  the  desired 
amount  of  flexibility  in  the  retrieval  language  as  well  as  the  ultimate  length 
of  the  fixed  records.  If  several  users  at  remote  sites  are  to  have  independent 
access  to  the  data  base,  there  must  be  terminals  at  these  sites  as  well.  Once 
the  hardware  availability  has  been  determined,  the  final  design  details  can  be 
drawn  - then  it  is  a matter  of  collecting  the  data  and  recording  it  on  the 
storage  medium.  To  be  useful  as  a working  tool,  the  data  management  and  retrieval 
system  must  then  be  maintained  in  a constant  state  of  timeliness,  and  there  must 
be  a continual  program  of  refinement  and  updating.  This  means  that  notice  of 
changes  to  the  system  be  automatically  sent  from  the  originators  of  such  changes 

to  the  data  base  librarian. 

Conclusions 

The  steps  towards  implementing  an  effective  Strategic  c3  data  management 

and  retrieval  system  have  been  outlined.  The  exact  nature  of  the  system  will  be 

dependent  upon  the  hardware  to  be  employed.  However,  even  a rudimentary  system 

'could  lead  to  improved  understanding  of  the  Cu  system.  If  used  to  advantage,  the 

data  base  will  not  be  a mere  catalog  of  recorded  information,  but  can  present 

the  user  with  an  analytical  tool  of  Strategic  C&  status  and  history.  This  tool 

can  also  be  used  by  researchers  at  locations  remote  from  the  host  computer. 

Extending  this  idea  further,  it  could  also  be  used  as  a standardized  base  of 

readily  available  inputs  for  all  Strategic  clt  modelling  and  simulation  efforts. 
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